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Green's and allied theorems : a historical sketch. 

By George A. Gibson, M.A. 

[Abstract.] 

The chief purpose of the paper was to indicate the rise of 

transformations of the type dxdydz = Vcosoc^S where the 

integral in the first member of the equation is taken throughout a 
closed surface, and that in the second member over the surface, a 
being the angle made with the axis of x by the normal to the 
element c?S drawn outwards. It is on this transformation the 
analytical proof of Green's theorem depends, and it was shown to 
have been employed in various forms by Poisson, Duhamel, Gauss, 
and others, before Green's essay was generally known on the 
Continent. It may be observed that the essay was published at 
Nottingham in 1828, and seems to have been unknown to con- 
tinental mathematicians till its reprint in Crelle's Journal, vols. 39 
(1850), 44 (1852), and 47 (1854). The following references were 
given in the paper: — 

Lagrange, in the Mdcanique Analytique (2nd edition, 1811), 
Part I., sect, vii., arts. 29, 30, gives the transformation : — 

^X'{8x'.dydz + 8y',dzdx + Bz'.dxdy) = Q A.'(cosa . Ex' + cos/3'. 8y' + 

cosy', 8z')d^, ds^ being an element of surface. 

Laplace, in the Supplement a la TMorie de V Action Capillaire 
(which forms a supplement to liv. x., Part II., of the Mdcanique 
C^ste, published 1806)— Oeuvres CompUtes, 1880, t4. pp. 428-432 
— transforms the integral 

r [dxdyi^ P j^ ^ ^ \ 

taken over the area of a section of a cylinder, whose generators are 
parallel to the axis of z, into the integral ± ^ Z — taken along 

Iv 

the boundary of the section, the + sign holding for the part of the 
curve convex to the axis of oj, the - sign for that concave to the 
same axis. 

Gauss gives a series of remarkable theorems, closely related to 
the transformation in question, in the introductory articles of his 



" Theoria Attractionis corporum sphseroidicorum ellipticorum homo- 
geneorum," Comment Soc, Beg, GotL, vol. ii. (1813), Werke, Bd. V. 
Great care is taken to make the proofs quite general. Some of 
these theorems will be found in Williamson's Integral Calculus 
(3rd edition), arts. 192, 193. Transformations of a like character 
occur in Gauss's Memoirs, "Principia Generalia Theoriae Figurse 
Fluidorum in Statu Aequilibrii," Comment Soc. Reg, Gott, vol. vii. 
(1830), Werke, Bd. V., and "Allgemeine Lehrsatze in Bez. auf 
die . . . Anziehungs-und Abstossungs-Krafte," ResuUate des Magn, 
VereinSy 1840, Werke, Bd. V. The proof of Poisson's equation, in 
§§9, 10 of the latter memoir, depends really on the transformation 
of a volume integral of the kind mentioned into the difference of a 
volume and of a surface integral, the manner of carrying out the 
transformation being quite special. 

Poisson in various memoirs uses with 'considerable effect the 
transformation employed (later) by Green, and establishes in his 
earlier memoirs particular cases of Green's theorem. In his 
"M^moire sur la Theorie du Magnetisme," MSm, de VInsty t.v., 
1826, read February 2, 1824, pp. 294-298, he proves the equation 

= r [k'ia'cod' + /J'cosm' + y'cosn')— 
J J P 

- r r [l d{'^'^') . d{P'k') d{y'k') \ d^'dy'dz' 
]\]\ dx' ^ dy' ^ dz' I p 

where p^ = (as - xy + (y- y'f + (» - a')' and the volume integrals are 
taken throughout a closed surface and the double integrals over the 
surface. He considers the possibility of a line cutting the surface 
in more points than two, but makes no reference to Gauss's memoir 
of 1813. In § 18 of the same memoir he discusses the case in which 
1//3 becomes infinite within the limits of integration. 

Another important example from Poisson is furnished by § 89 
of the sixth chapter of his Thdorie de la Chaleur (Paris, 1835 ; a 
date later than that of the publication of Green's essay and of the 
memoir of Duhamel referred to below). He arrives at the following 
equation : — 



m 



cuVdxdydz 



= PA;( -— cosa + -r-cosiS + — cosylc^S 

J J \dx ay dz / 

- uk\ — cosa + — cosp + — cosy laS 

]] \dx dy ^ dz ^ f 



+ 

The last equation is in form the extension of Green's theorem given 
in Thomson and Tait's Natural Philosophy, vol. i., p. 168. 

Another interesting example from Poisson occurs in a " M^moire 
sur I'Equilibre des Fluides," M^m, de Vlnst,, t.ix, 1830, read Nov. 
24, 1828, pp. 22-24. 

Ostrogradsky, in a memoir with the title " Note sur la Th^orie 
de la Chaleur," M^7n. de VAcod. de St Fetersbourg, 6s. t.i (1831), 
read Nov. 5, 1828, p. 129, establishes the theorem 

r/|? + J+^k= r/PcosA + Qcos/x + Rcosi/)s 

where w is a volume-element, s a surface-element, and P, Q, R the 
values of JO, q, r at s, 

BuhameJy in Note 1. appended to a memoir on the "Theory of 
Heat," published in the Journ, Ecol, Polyt, t.xiv, cah. 22 (1833), 
gives, at pp. 67-71, a proof of the equation 



+ -^r-=r + —^-ir-]dxdydz 



J \da^ df dz" f 

-h UlZ-f-- -Hm— — -hn— -— IctS 
J \ dx dy dz f 

-H 



dy 

I-—- + m + n laS 

dx dy dz / 



Z, m, w, being the direction cosines of the outward normal at c?S. 

Lamdy in the same volume, p. 204 and p. 231, gives investiga- 
tions of a similar character. 

Sir W, Thomson, in Article XII. of his "Electrostatics and 
Magnetism " (reprinted from the Camh. Math. Jour., Nov., 1842, 
and Feb. 1843), furnishes examples of the analysis with which the 
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paper dealt, and which, he says, was suggested to him by the 
analysis used by Poisson in the article of his Thdorie de la Chalen/r, 
quoted above. 

Questions of priority are usually somewhat difficult to answer ; 
but while it seems clear that the theorem generally quoted as 
Green's was given independently . of Green, yet the importance 
which he rightly attached to it, and the splendid use to which he 
put it, amply justify us in keeping to the customary mode of 
citation. 



Some now Properties of the Triangle. 
By J. S. Mackay, M.A., LL.D. 

[The substance of this paper will be included in Dr Mackay's 
paper on the triangle in the first volume of the Proceedings, now 
about to be published.] 



Second Meeting, 13th December 1889. 



R. E Allardice, Esq., M.A., Vice-President, in the Chair. 



A special case of three-bar motion. 
By Professor Steggall. 

The questions involved in the consideration of three-bar motion 
have attracted a good deal of attention (Proceedings of Mathematical 
Society of London passim, and elsewhere) ; but I am not aware 
of any complete account of the figures that can be derived from such 
a motion. The present paper gives a complete list of all the different 
kinds of curve that are obtained by a tracing point at the middle of 
the middle bar, the two outer bars being equal. 

It may be advisable to briefly obtain the general equation to the 
curve traced by any point on the middle bar, without any condition 
of equality in the lengths of the other two. 

Let 2a be the distance of the fixed centres, 6, 2c, d the lengths 
of the three bars in order, h the distance of the tracing point from 
the middle of the middle bar measured from the bar 6, 6, <^, xp the 



angles which the bars in order make with the line joining the fixed 
centres : take this line as axis of x, and a perpendicular through its 
middle point as that of y. 

We have at once, on reference to a diagram, 

6cos 6 + (h + c)cos ff> — a-\-x^ 
6sin ^ + (^ + c)sin <^ = y, 
(fees ^ + (7i - c)cos <l> = x-a, 
dsin \j/ + {h — c)sin <^ = y. 
Whence 

- 2(h + c){(x + a)cos <^ + ysin <^}, 

(P=f+{x-ay+{h-cy 

- 2(^ - c){{x - a)cos <^ + ysin <^}. 
Whence 

i{(hx + ac)co8 <t> + hysia<{>} = 2{a^ + f + a:' + h^ + (^)-b^-<P, 
4 { (ex + ah)cos <f> + cysin<^} = iax + ihc -b^ + cP. 
Now the eliminant of 

Pcos <^ + Qsin <f> = 'R 
P'cos <^ + Q'sin <^ = R' 
is 

R^F^ + Q'2) + R'2(p + Q2) - 2RR'(PF + QQ') = (PQ' - QP^. 
Whence, on substitution, we easily obtain 
(calling 2a2 + 2h^ + 2c'-b^-(P, A, and Uc + (P-b\ B) 

r«. 4c« 

- r*x . Sack 

+ r" . 4{c(cA - hB) + aVv"} 

-r^Q^ . 16aV 

-r'x. 4a{B((r*-A2) + 4a2cA} 

+ aj^ . 32a'cA 

+ t" . {4a2A(AA - cB) + (cA - 7aB)« - IGa^c* - A*)^} 

+ ar^ . 8a2{2aV - c(cA - KB) - A(/iA - cB) + 2(c' - A^} 

+ X . 2a(AA - cB)(cA - KB - 4a2c) 

+ a2(AA - cB)« = 0. 

It may be worth while to write this at full length : the result is 
r«. 4c'» 

- r^x , Sack 

+ r* . 4:{a\Jv' + 2c') + 2(^(0" - h') - b^c - h) - (^^0(0 + A)} 

-r'a^. 16aV 

+ r^aj . 4a{(A2 - c2)(4Ac + c^=^ - 5^) - 4a2cA} 

+ a^ . 32a'cA 



+ [2c(c2 - A^) - b\c - ;*) - (Tic + h)f} 
+ aP . 8a%b\c - hy + (P{c + hy - 2a2A«} 
-x,2a {2h(a^ + h^ - c") + h\c - A) - c^(c + h)) 
{2c(a» + A^ - c^) + 5'(c - A) + (^^(c + A)} 
+ a2{2A(a2 + A'' - c) + 5«(c - /i) - (^(c + A)}' = 0. 

The same equation has been obtained by a slightly different 
method of elimination by Professor W. W. Johnson {Messenger of 
McUhematicSy vol. V., 1875, p. 50). 

The equation is so unmanageable in this form that some special 
assumption seems necessary in order to be able to trace the sextic 
curve ; and the assumption in this paper is that b = d, h = ; in 
other words, the tracing point bisects the free link, while the fixed 
links are of equal length. 
The curve now becomes 

r« + 27^(a» + c^ - 62) - ^r'x'a' 
+ r'iaXaJ' - 2c* - 26«) + (c* - by} 
+ Wb^^ = 0. 
This equation may be readily solved to give r in terms of ^, the 
vectorial angle, and the result (which may be easily obtained in 
other ways) is 

r» = a2cos2 ^ + 6^ - c^ + 2asin Jc" - a^cos^ 6, 
or sin (9 = (a^ - c* + 6* - r')l2a J{b^ - r*). 

The minimum values of r- are when cos = 0,sin 6 = l,cos2 ^ = - 1, 

and in this case T^ = b^-{a + cy. This expression being only positive 

if b>(a + c) the division is at once suggested into the cases 

I. b>a-\-c 
II. b = a + c 
III. b<a + c. 
Before proceeding to the separate cases, we may notice that the 
maximum and minimum values of r^ are given by the equation 

rdr/dO = = - a-sin2 6 + acos e{c' - a'co&2 6)/ J{c' - a^cos- 6), 
whose only solution is cos ^ = ; and those of y^ by the equation 
sin ^cos ^{4a2cos2 d^-b^-c"- 3a* 
± asin ^(3c* + a* - 4a-cos- 6)1 J{c^ - a*cos« ^) } = 0, 
whose solutions are cos ^ = 0,sin ^ = 0, and 
Sa*b^cos* + aW{b^ - lOc^ - 6a-}cos* 6^ {(a - cy+h''a} {{a+cy - Wc] = 

... ... ... ... ... ... IXJ 

or 16a*cos* 6 = 10c* + 6a* - 6* ± {5* + 2(c* - a*)} ^1 + 8(c* - a*)/6* (2) 



I. b>a + c. 

We notice that if a'^c, can have every value, otherwise the 
maximum value of cos is c/a : we therefore divide this case into 

(i.) a<c, (ii.) a = c, (iii) a>c, 
(i.) a<c. 

Here may have every value giving always two values of r* 
except at ^ = 0. We get two detached but intersecting ovals which 
must be separately described, as may be at once seen on drawing 
a diagram. 

On referring to equation (1) we see that if h^>{a + cy/c, the pro- 
duct of the values of cos' is negative, and the negative value of cos' 
is inadmissible, while from equation (2) it is clear that the positive 
value of cos' is greater than unity. Thus the only maximum and 
minimum heights are on the axis of y. The form of one branch of 
the curves is given in figure i. In this, and in all the other figures, 
the three numbers affixed denote the values of the constants a, 6. c, 
taken in that order. As a rule the standard value of 6 is 12, but 
for clearness it is taken of various convenient magnitudes. 

If 5'<(a + c)'/c, both values of cos'^ are positive, and the smaller 
value is less than unity. This gives two maximum values of y on 
each loop, symmetrically situated with regard to the axis of y; 
while if 6' = (a + cy/c these two maxima coalesce with the minimum 
between them, giving a flattish figure to the curve (see figures ii. 
and iii.), and a point where the tangent meets it in four coincident 
points. 

(ii.) c=ia. 

In this case, since the two loops in general make (on opposite 
sides) an angle with the axis of a;, where they cut, of 

tan-^{(a' + 6' - c^^a J^c" - a')} 
the two loops cut the axis of x at right angles, and therefore 
touch and may each be continued from the other, since the lines 
BO, AD when they lie along AB may be made to cross or not cross 
in continuing their motion, as may be desired. 
The equation reduces to 

r* = 6' - 2a'sin' ± 2a'sin J{\- cos' 0\ 
which, by taking the top half of each with the bottom of the other, 
becomes 

r' = 6' 

r* = 6' - 4a'sin' 0. 



The latter is the inverse of an ellipse with regard to its centre, 
since 6 > 2a ; and it becomes two circles if 5 = 2a. 

This curve being written 

y2={54_(52_2r2n/16a2 
the least value of r^ is 6^-4a-,^ = 7r/2; and if this is >b^l2y the 
greatest value of 2/* is when r^ has this least value ; this requires b^ 
to be >8a' ; but if 5'<8a^, and therefore b^/2>b^ - ia\ r^ can be as 
small as 672, and this then gives a TtwoAinwm, value of y — viz., 
b^l^a^sin^O = 678a* ; and a minimum value of y when = 7r/2. Also 
when 6' = 8a' these two coalesce giving a masked point of inflexion 
on the axis of y. See figures iv., v., vi. 

It is interesting to notice that as a approaches c two disconnected 
equal portions of the loops become the circle, and the inverse of the 
conic ; so that each loop degenerates into a semi-circle, and half of 
the inverse of the conic. 

(iii.) a>c. 

In this case cos is limited, its maximum value being c/a. 

If b^>{a-\-cfjCy we find that the values of cos'* in equation (1) 
are one positive and less than d^ja^ (as may be at once verified by 
substituting cja for cos^ in the expression in the left which is then 
positive, while it is negative if cos is put equal to zero), and the 
other negative. This gives us one minimum value of y\ and its two 
positions are on the lower part of the upper loop and vice versa. 
See figure vii. 

It is noteworthy that to pass from the case a = c^ to the case 
a>c, we may regard the upper half of the circle as combined with 
the upper half of the inverse curve, whereas in passing to the case 
a<c, as seen above, the upper half of the circle is combined with 
the lower half of the inverse curve. 

If 5* = (a + cYJc the equation for cos- becomes 

8a«cos ^e + (6= - 10c' - 6a2)cos '^ = 0, 
giving cos '^ = or (a' + 3c')(3c - a)/8a'c, 

the last value is always less than <?la^y and we thus have, in this 
case, an additional zero value of cos' 6, and the tangent at the vertex 
meets the curve in four coincident points. If 3c>a there are two 
other symmetrical minima on each loop (fig. viii.), if 3c = a these 
also move up to the vertex where the tangent now meets the curve 
at six coincident points (fig. ix.), and if 3c>a, these minima dis- 
appear {^g. X.). 
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The equation to the curve of £g, ix. is 
iB« + 3aj*(y2 - 48c*) + Sa^if - iSc^)(f - 36c«) + y'if - 48c»)(y» - 60c*) = 
from which it appears that the lines y' = ^Sc* each meet^the curve in 
six coincident points, and the lower part of the curve approaches very 
closely to a straight line. 

If b^<{a + cy/c clearly no other maxima or minima exist unless 
6*< lOc* + 6a' (fig. xi.) ; if 6*< 10c* + 6a*, we have to divide the cases 
into, from equation (2), (i.) 6*>8(a* - c*) in which case both values if 
cos* are real and admissible giving four )symmetrical) maxima or 
minima heights on each loop (fig. xii.); (ii.) 6* = 8(a* - c*) in which case 
these coalesce two and two giving a point of inflexion with a hori- 
zontal tangent (fig. xiii.); (iii.) 6*<8(a*-c*), in which case the 
points of inflexion are left but the horizontal tangent becomes 
oblique (fig. xiv.). 

II. b = a + c. 

This case naturally divides like the last into 

(i.) a<c, (ii.) a = c, (iii.) a>c; 
and in each case 



r* = 2(a*cos ^6 + oc ±asin 6 Jc^ - a*cos *e). 




or r = Jail + cos 0){c + acos 6) ± Ja{l - cos 0)(c - ocos 6), 

The equation to obtain the maximum height is now, besides 

cose = 0, 1 6a*cos* ^ = 9c* - 2ac + 5a* ± (3c - a) J{c + a)(9c - 7a). 

Case (i.) a<c. 

Here the upper sign gives a value of cos* $ greater than unity, 
and the lower sign gives a positive root less than unity, as may be 
seen from the more general case I. (i.), or by comparing (9c* - 2ac + 
5a*)*, and (9c* - 2ac - 11a*)* with (3c - a)*(9c* + 2ac - 7a*). 

Hence in this case there is one maximum height besides that on 
the axis (fig. xv.). The two minima heights on the axis here become 
zero, the two loops meet one another at the origin with a common 
vertical tangent there, so that there is a choice of path at the central 
point. 

Case (ii.) a = c. 

Here the pairs of tangents at the points of intersection on the 
axis of Xf which were inclined to that axis, become coincident, and 
we simply get three circles (fig. xvi.). 
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Case (iii.) a>c, 

In this case we have, as in I. (iii.), the condition cos ^O'^^t^ja^ ; 
subject to this both values of r^ are real. For the maximum height 
we have the same equation as in case (i.), namely, 16a'cos^^ = 9c'- 

2ac + ba^ ± (3c - a) J{c + a){dc - 7a), and the values of cos ^0 are real 
if 9c is equal to or greater than 7a, or if 3c = a. Now if 7a<9c, it 
is easily seen that both values of cos ^0 are admissible, and we have 
(fig. xvii.) four (symmetrical) maxima and minima : if 7a = 9c, these 
coincide pair and pair, giving a point of inflexion with a horizontal 
tangent (fig. xviii.). 

If 7a>9c (fig. xix.) these points disappear, while if a = 3c, the 
values of cos*^ though real give imaginary values of r*, and the 
figure resembles that just referred to. 

III. b<a + c. 
Writing our original equations in the form 



r' = aJ' + b^-c'^ 2asin ^(osin 0± J(^- a^cos ""0), 

sin e = {c2 - (a^ + 6^) + r2}/2a Jb^ - r^, 
we shall find a sub-division c^> = <a* + 6* convenient ; but we shall 
keep to the other division of a > = < c as before. 

Case (i.) a<c. 

If a^ + b^<c^y sin increases with r, and lies between {c^ — a^ - 
b^}/2ab and 1, and each value of sin gives only one positive 
value of r* : the least value of r^ is zero, sin ^ = {c* - a^ - b'}/2ab, and 
the greatest is 6^ - (a - c)% sin ^ = 1 : the curve consists of a simple 
figure of 8 (fig. xx.). 

If a* + 6^ = c^, the two loops touch with the axis of x for common 
tangent meeting the curve in six coincident points as in tig. xxi. : the 
equation in fact reduces to 4ay (6^* - r^) = r^ = (or* + q/^y. 

If a^ + 6^>c^, the values of r^ are sometimes both positive,, some- 
times not : the limiting value of is clearly given by 

a2 + 62 - c2<^ 2asin ^(asin + Jc" - a^cos ^6), 
or, in other words, sin must not be greater than the value of that 
makes r* vanish : thus one available value of r* exists for all values 
of 6 and two values form sin ^ = to sin ^ = (a'* + 6^ - c^)/2ab. 

The maximum height is found as before, and refers to the loop 
(fig. xxi). In this case the curve is described by a single operation 
without any choice of direction at any point, or break of continuity. 
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Case (ii.) a = c. 

Here we get at once the equations 

r^ = h\ and r^ = 6* - 4a'sin *^, 
which last, since 6 < 2a, is the inverse of a hyperbola that becomes 
rectangular when 6^ = 2a'* (fig. xxiii.). 

Case (iii.) a>c. 



In this case cos 0<c/a,sin 6> J\- c^/a\ and since we have 

sin e = {(a^ - c^) + (6^ - r')}/^CL Jb^ - r«, 
we see that sin decreases with b^-r^ until, if possible, 6^ - r* = 
a^ - c\ at which point it begins to increase again until sin ^=1. 
When r = 0,sin 6 = (a^ + b^- c^)l2ab, which is always an admissible 

value : and between sin^ = (a^ + b^- c*)/2a5and sin ^ = ^1 - c^/a there 
are always two values of r^, it only remains to see whether r^ can 
ever become 6^ 4- c* - ^* : this of course depends on whether a^< = > 
b^ + c\ 

If a^<{b^ + c^) there are two real values of r between these limits, 
and one beyond, t.e., from sin ^ = (a^ + 6^ - c^)/2ab to sin ^ = 1. On 
referring to the equations (1) and (2) we find that if 6'>8(a'-c') 
(and therefore b^ + c^>a^) there are four symmetrical maximum and 
minimum heights (fig. xxiv.); if b^=^S{a^-c^) there are only two, 
through the coincidence of two pairs forming a point of inflexion 
with a horizontal tangent (fig. xxv.), while if 6^<8(a*-c^) there are 
none (fig. xxvi.). 

If a' = 6^ + c^, the two limits coincide and we have a curve with a 
very approximately straight portion near the node : this, in fact, (fig* 
xxvii) is Watts' parallel motion. 

If a* > 6* + c\ the value of 6 decreases till r^ = 0, but cannot decrease 
further (fig. xxviii.) through the disappearance of the limit sin ^ = 

In all the descriptions of these curves the upper half only has 
been referred to ; and from the obvious symmetry it has only been 
necessary to discuss values of less than 7r/2. It is clear that we might 
from the beginning have proceeded under the heads a<c,a = c^a>c; 
but some consideration has led to the adoption of the process of this 
paper. There being two independent ratios, there are numerous 
cross connections between the curves, but these are best seen by the 
help of the diagrams. 
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On some theorems in the theory of numbers. 
By B. E. Allardige, M.A. 

The number of groups of n which may be selected from 2n is 
27i(2n - l),..(w + l)/n ! But make the 2n into two groups of n, and 
select r out of the first and n-r out of the second. This gives 
[w(n-l)...(n-r + l)/r!] + [w(n-l)...(r+l)/(n~r)!] ways of thus 
making a group of n. Hence 

♦2n(2n-l)...(n + l)/n! = l+w« + [n(n-l)/2!]«+ (1). 

.-. 2n(2n-l)...(n+l)/n!-2 

= n«{P + [(n-l)/2!p + [(n-l)(n-2)/3!]»+...} 
= 7i^{P,^ + P,«+...+P„?.J(say). 

We shall now show that Pj* + P2* H h P„Li is divisible by w, if 

n be prime. 

Py = P, (mod.n.), if r + 8 = n, but not otherwise. For if 
P^ = P, (r>«), then 

(n-l)(n-2)...(n-r+l)/r!-(/i~l)(/i-2)...(n-« + l)/«! = 0; 
(n-l)(n-2)...(n- g+l) / (n-g)(n~g-l)...(n-r+l) ,\^Q, 
•* ^! I («+!)(« + 2)... i^ /~ ^ 

.-. (n-s)(n-«+i)(/i-« + 2)...(n-r^)-(s+l)(ff + 2)...r = 0; 
.-. ±»(« + l)(« + 2)...(r-l)-(« + l)(s + 2)...r = 0; 
.-. -(« + l)(« + 2)...(r-l)( + s + r) = 0; 

and this is true if r + « = w (otherwise obvious) and not in any other 
case. [If r + « = w, then r - « = w - 2«, which is odd, and the lower 
sign is to be taken where the double sign is printed.] 

It is obvious that P^ + P, is not divisible by n ; and hence if we 
divide Pi, P2,...P(«_i)A! ^7 ^> we must get for remainders either 1 or 
n - 1 and either 2 or n - 2 and so on. 

Now since (n - rf = n^- 2nr + 7^ = 7^ (mod. w), we must have 

P,« + P22+...P(Li)/2=l' + 3« + 52...+(7i-2)« 

= »(«-l)(n-2)/6; 

which is divisible by w if w be any prime except 2 or 3. 
From this, and the identity (1), it follows that 

(2n-l)(2w-2)...(7i+l)-(n-l)(n-2)...l = 0(mod. »'). 



* The nee of this identity was suggested to me by Professor Tait. 
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We shall next show that 11+ — + — +... + -i(n-l)! is 



044-^"^n--T)("-l)' u 



divisible by n^ * 
We have 

\ 2 3 n-\r ' \l.(n-l) 2.(n-2) P ' 

Hence we have to show that I—- __ + __ tt- + ...|(w~ 1)! is 

\l.(w-l) 2.(w-2) r ' 

exactly divisible by n. 

Assume (n - l)!/l.(w - 1) = Oj, (w - l)!/2.(w - 2) =: C4, &c. 
Then 



(r + l)V+i - ^^r = (*• + l)*(w - 1 i)/(»' + l)(n - r - 1) - r»(w - 1 !)/r(n - r) 

= {(r + l)/(n - r - 1) - r/(w - r)}(w - 1)! 
= (w!)/(n - T)(n - r - 1) = (mod. n). 
/. (r+ 1)^+1 = r2a^ 

^(r-l)V-i 



I • • I 



= 1 (by Wilson's theorem). 

Hence we may write 

"1 = ^f4 + 1 
2*03 = n/ig + 1 



••• ••• ••• 



m'a^ = w/i^ + 1 (where m = (n - 1 )/2) 

. •. (mKflsa, = P.w + (m!)«( l/P + 1/2^ + . . . + l/m^). 

Now, if we assume (m !)/r = a„ we may easily show that a^ + a, is 
not divisible by n, and hence that 

ai' + a2^+... +aj, = (mod. w), 
which proves the theorem. 

Consider again the result 

A = (2n-l)(2n-2)...(n+l)-(n-l)(n-2)...l=0(mod. »'). 
This gives 

A = (w + w-l)(7i + w-2)...(n+l)-(w~l)(n-2)...l 

where jp„_3 = (n - 1 !)2( 1 lr8){^ =N «) 

/>„-» = (^^l!)2(l/r). 



* Compare a paper by Mr Leudesdorf, in the Proceedings of the Lond, 
Math, Soe. for 1889, p. 199 — a paper which I did not see till after the above 
was written. 
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No^ Pn-Jt ifi divisible by n', and hence p„^ is divisible by n. 

This theorem may also be proved in the following manner:- 



We have 2(n - 1!)2(1M 



= {(w- i!)/1.2 + (n-l!)/1.3 + ...+(n-^l!)/l.(n-l)}( = P0 
+ {(^ir^\)l2.l + (wrn:!)/2.3 + ... + (irri!)/2.(n - 1)}( = P,) 

iT ••• ••• •*. ••• 

Now consider the terms of P„ namely, 

{^l \)lr.l,{^l !)/r.2, (w^l \)/r.{r - l),(^r^l !)/r.(r + 1), &c. 

No two of these can be congruent ; and 

(w - 1 l)/r.p + (n - 1 !)r.(n - j^) = nljr.p.^n -p) = (mod. n). 
Hence if we divide each of the terms of P^ by «, we get as 
remainders all the numbers 1, 2, 3...n-l, with the exception of 
that number which is complimentary to a^ where a^ 

= (w - 1 \)lr.(n - r) (mod. n). 
Hence the sum of all the remainders in 22(n - 1 !)/r« 

= (n- 1)(1 + 2 + ... +n^) - 2(01 + 03+ ... +a(„_„) 
= (n-l)'n/2-2{l/l.(w-l) + l/2.(w-2)+...} 
which is divisible by n. 

The theorem that the sum of the reciprocals of the numbers 

1, 2,...w- 1, is divisible by n*, when n is a prime, may be extended 
to the sum of the m** powers of these numbers, where m is an 
integer, positive or negative. 

Let S^=l"* + 2"*+ ...(w-1)"*; it being understood that if m is 
negative ( = - Q, the sum of the powers is to be multiplied by 

(n - 1 !)', so that it may be made integral. 
Since, when n is prime the equation 

(a; - 1 )(a; - 2) ... (a; - w-1 ) - a*^^ + 1 = (mod. n) 
has (w-1) incongruent solutions, each co-efficient is divisible by n. 
Hence, if m is positive, S^ is divisible by n, unless m is a multiple 
of (w-l). 

Suppose now that m is an odd positive integer and w=*= 2 ; then 
2 S^ = 22a"» = 2(a"» + n-o"*) 

= 2{ a"» + n"* - rrfljnr-^a + . . . + ^Cino*^^ - a"*} 

and S^_i = (mod. n), unless m - 1 is a multiple of n - 1 ; 

S^ = (mod. n^), unless m - 1 is a multiple of w - 1 ; 

and the theorem is true even in this last case if m is a multiple 
of w. 
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Now consider 



S^^={l + l/2- + l/3"'+... + l/(n-l)-}(n-l!)"'. 
We have 



»-i 



2 S_^ = (n - 1 !)'»2{ l/r"* + l/(n - r)~} 

— (n - 1 i ) ^ ^/ rzr- 

-(n-l!)2 -;;.(;^::^:y;^ 

We have thus to show that ^Ci(w - l!)"*2{l/r(w-r)"*} is divisible 
by w. We shall suppose, for the sake of clearness, that m is less 
than n ; but the following method will be applicable, even if in be 
greater than n. 

Assume (w - 1 !)"*/{ (»i - r)r"*} = a^ (mod. n) 

(n-1 !)"* = a^(n - r)r*^ 

a^(n-r)r^=-l (by Wilson's theorem). 
Now since (n - rY~^ = r^~^ = 1, we get 

a^ = - r"-"*"^(n - r)"~^ 
a^ = r»^»t-ir»-2 (mod. n) 

Hence 2a^ = 27'**-^2 

= 0, if w-m-2 + 0. 

It follows that S_^ is divisible by Wj, m being subject to the 
restriction n- 171 — 2 be not zero. If we remove the condition that m 
is to be less than w, we shall easily find that the general restriction 
as to the value of m, is that m + l must not be a multiple of n - 1. 

In the paper referred to before in a footnote, Mr Leudesdorf 
considers the case where n is not prime and S^ denotes the sum of 
the m'* powers of the numbers less than n and prime to it. His 
method however cannot be considered rigorous, as it involves the 
use of divergent series. 



Note on normals to conies. 
By R. H. PiNKERTON, M.A. 

1. The following condition may be new ; it does not appear in 
any of the books : — 

The condition that the straight line 
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lx + my + n=^0 (1) 

should be a normal to the general conic, 

8 = aa^ + 2hxy'\-by^ + 2gx + 2/y + c = (2), 

referred to rectangular axes, is 

^aP + 2hlm + brri')^^{P + my (a), 

where 2 is written for 

A^^ + Bm" + CnJ" + 2¥mn + 2Gnl + 2mm, 
and A, A, B, C, F, G, H have their usual meanings. 

We know* that the equation to the tangents to S at the points 
where (1) cuts S, is 

S2 = A{lx + my + n)\ 
or a?{a2 - N^) + 2xy{h:2 - Mm) + f(h2 - Am*) 

+ terms of lower degree = (3). 

Now the straight line (1) will be normal to S if it is perpendicular 
to one of the straight lines (3). The condition for this is 
Z»(a2 - M^) + 2Zm(A2 - Mm) + m\h2 - Air,*) = 0, 
or ^{al^ + 2hlm + hm^) = A(Z* + 2ZW + m% which is (a). 

2. If the line (1) be given as passing through a given point 
(a, /3) its equation will be 

Ix + my - (la + m/B) = 0. 
If in the condition (a) we substitute -la- mP for w, we shall get 
the following biquadratic in l/m — 
(a^ + 2hlm + bm^) x 

[AP + Bm* + C{la + mpy - 2Fm{la + m/3) - 2Gl{la + mfi) + 2mm] 
=:A{P + my, 
or l'[a{A + Oa^ - 2Ga) - A] 

+ 2Pm[a{Cal3 - Fa - G^ + H) + A(A + Ca* - 2Ga)] 
+ Pm\a{B + 0^- 2F/3) + b{A + Ca^ - 2Ga) 

+ a(Ca^ - Fa - GjS + H) - 2A] 
+ 2lm\b(Cap - Fa - Gj8 + H) + h{B + 0^8* - 2¥P)] 

+ m*[5(B + C/?2-2F^)-A] = (/]?). 

This biquadratic gives the directions of the four normals which 
can be drawn from the point (a, p) to the conic S. 

3. If (», y) be any point on any one of the four normals from 
(a, P) to the conic S, then 

l(x - a) + m(y - /5) = 



♦ Salmon's Conies —Fitth Edition— Art. 374. 
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if Ijm be one of the roots of (/?). 

Hence ^/{y~P)= -m/(x-a). 

If then for I and m we substitute in (P) y - P and - (a; - a) re" 
spectively, we shall obtain the equation to the four normals which 
can be drawn from {a,P) to the conic S. 

4. It may be worth considering what the condition (a) becomes 
for the circle — 

ar* + y2 + 2^a; + 2/y + c = 0. 
Here a = 5 = l;A = 

A = c-/^,B = c-(^,C=l, 
F=-/, G=-^, H=/^. 
Hence in this case (a) becomes 

(P + 7n%P{c "f) + m\c - ^2) + ^2 . 2/mn - 2gnl + 2fglm] 

Dividing by P + m\ we obtain 

c{P + m^) -J^P - g^rn^ + w^ - 2/mn - 2gnl + 2fglm 
= c{P + w?) -fl^ -frn^ "fJ^- g'mS 
or /Z^ +J^m^ + 2^Zm - 2n(^Z -{-fm) + w' = 0, 

or (^Z +/m - nf = 0. 

Hence the condition that Ix + my + w = should be a normal to the 
circle is gl +/m -n = 0» 

But this is the tangential equation to the point {-g, -/), the centre 
of the circle. 



Third Meeting , 10th January 1890. 



George A. Gibson, Esq., M.A., Ex-President, in the Chair. 



Note on a carioas operational Theorem. 
By Professor Tait. 

The idea in the following note is evidently capable of very wide 
development, but it can be made clear by a very simple example. 
Whatever be the vectors a, /?, y, 8, we have always 

V. Va^VyS = aS.^SyS - ^SS.ayS. 
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Bat vector operators are to be treated in all respects like vectors, 
provided each be always kept hefort its subject. 

Let o- = t^+yiy + *i 

where ^,i7,{ are fanctions of x^y^z ; and let 

.d . .d ^ d 

as usnaL Also let o-^, y^ be their values when x^, y^, % are put for 

Then by the first equation, attending to the rule for the place of 
an operator, 

V. Vyo-VyiCTi = yS-oyjO-j - S{yi(riy)a-. 

If we suppose the operations to be completed, and then make 
a^ = ae,yi = y, 2^ = 2, the left-hand member must obviously vanish. So 
therefore must the right. 

That is : — yS-oyiO-j = S(yicriy)a- ; 

if when the operations are complete, we put o-i = o-,yi = y. 

In Cartesian co-ordinates this is equivalent to three equations, 
of the same type. I write only one, viz. : — 



d^ 
dx 



^ 1 f 

d^ d_ d_ 

dxi dffi dzj 

f 1 Vi Ci 



d^ 
dxi 
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dz 
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if, afier operating, we put a^ = *i f i = ^> <fcC'i ^^ 



On a property of odd and even polygons. 
By R E. Allardige, M.A. 

The property referred to comes to light on consideration of the 
problem, "To inscribe in a given n-gon the n-gon of mininnimri 
perimeter.'' 



Tbiangle. 

Let us consider first the case of the triangle (fig. 29). If ABO 

is the given triangle and DEF the inscribed triangle of minimum 

perimeter, it is obvious that we must have l FDB = z. EDC( = a, say), 

L DEC = L FE A( = p), L EFA = l DFB( = y). This condition is 

satisfied if D,E,F, are the feet of the perpendiculars from the oppo- 
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site vertices ; and it is usually assumed that the triangle so obtained, 
the pedal triangle, must therefore be the triangle required ; in other 
words, the assumption is made that only one inscribed triangle may 
be constructed, whose sides shall be equally inclined to the sides of 
ABC. This assumption is as a matter of fact correct in the case of 
the triangle ; but it can hardly be considered not to require proof, 
as the corresponding assumption would be wrong in the case of the 
quadrilateral. 

We might proceed to determine DEF as follows : — 
We have a + )8 + C = tt, 

)8 + 7 + A = 7r, 

7 + a+B=7r; 
from which we get a = A, )8 = B, y = C. 

To calculate CD, put CD = a;, AE = y, BF = », then 

scsina — {h- 2/)sin)8, 

ysin^S = (c - z)^m.y, 

z^iny = (a - a:)sina ; 
whence asina = hsmfi - csiny + (a - a;)sina ; 

from which we get cc = (a^ + 6^ - c*)/2a ; which shows that D is the 
foot of the perpendicular from A on BC. 

Quadrilateral. 

Consider next the case of the quadrilateral (fig. 30), 
Let ^APS= ^BPQ = a; z.BQP= z.CQR = /?; etc. 

If now we try to determine a as we did in the case of the triangle, 
we get the equations 

a+)8 + B = 7r, /? + 7 + C = ^, 

y + 8 + D = 7r, 8 + a + A = 7r, 

whence B-C + D-A = 0; 

and the problem is therefore impossible unless the given quadrilateral 
be cyclic ; that is to say, it is in general impossible to inscribe in a 
given quadrilateral a quadrilateral each pair of consecutive sides of 
which shall be equally inclined to the side of the given quadrilateral 
in which they meet. 

If the given quadrilateral be cyclic there are an infinite number 
of inscribed quadrilaterals satisfying the given condition. 

The investigation that determined the position of the vertices in 
the case of the triangle determines the inclination of the sides in 
the case of the quadrilateral. 
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Using a corresponding notation in this case, we have as before 

a:sina = (b- 2/)sin/?, zsiny ~{d- to)sin5, 

ysinP = {c- 2;)sin7, twsinS = (a-x) sinS ; 

whence asina - bsiiiP + csiny - dsinS = 0. 

Let now a denote the angle subtended at the circumference of 
the circumcircle by the side a, jS' tte angle subtended by the side 
b; and so on. Then a = 2Rsina', and therefore 

sina sina - sinS'sinS = sin)8'sin)8 - siny'siny. 

But sina' = sin( A + 8') ; sin)8' = sin(0 + y'), 

sinS = sin( A + a), siny =sin(C + /?); 
. •. sin(A + 8')sina - sin8'sin(A + a) = sin(C + y)sin)8 - siny'sin(C + /?); 

. •. sinAsin(a - 8') = sinOsin()8 - y') ; 
and sinA = sinO, . *. sin(a - 8') = sin(/? - y'). 
Hence either a - 8' = tt - ()8 - y'), and therefore a + /? = 7r + y' + 6', 
which is impossible ; or a- 8' = 13- y' and . *. a - /? = 5' - y' ; 
and we have also a + )8 = tt - y' - 8', 

.-. a = 7r/2-y', /? = ^/2-5'. 

A geometrical investigation may also be given, in the following 
manner. 

Let ABCD (fig. 31) be the given quadrilateral. Take U any 
point in DA ; K the image of U in AB ; L the image of K in BC ; 
M the image of L in CD. Join MU, XL, WK, VU. Then it may 
easily be shown, that if ABCD is cyclic, XJVWX has its adjacent 
side equally inclined to the sides of ABCD. 

If U move along DA, the locus of K will be AR, the locus of L 
will be RS, the locus of M will be SM. 

Now z.SMU= z.SLX= ^VKA= ^VUA= ^XUD; .-. SM 
is parallel to DA. 

Now let U' be another point in DA; K', L', M', the points 
corresponding to K, L, M. Then UU' = KK' = LL' = MM'. Hence 
MM'U'U is a parallelogram ; the direction of UM is invariable ; 
and UM which is equal to the perimeter of the inscribed quadri- 
lateral is of constant length. [It should be noted that if the quadri- 
lateral UVWX be crossed, one of its sides has to be considered 
negative.] 

Having shown that the directions of the sides of the inscribed 
quadrilateral are invariable, and that the primeter is constant, we 
may make use of a particular case to determine the directions of the 
sides and the length of the perimeter. 
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Let one of the vertices of the inscribed quadrilateral (fig. 32) 
coincide with D, so that PQ is equal to the primeter. 

The triangles KCD and HAD are similar, and therefore the 
triangles KDH and CDA are similar ; hence 

KH:CA = KD:DC = sinO 
.-. KH = CAsinO 
. •. PQ = 2KH = 20AsinC = AA;/R, 
where h and k are the diagonals of the quadrilateral and R is the 
radius of the circumcircle. 

Denoting DMH by a, as before, we get 

7r/2-a= ^Q= z.KHD= Z.CAD; 
.'. a = 7r/2-CAD; 
that is, a is the complement of the angle subtended at the circum- 
ference by the side of ABOD opposite to that in which the vertex 
of a lies. 

It might be thought that there must in every case be some 
minimum inscribed quadrilateral. But if the given quadrilateral is 
not cyclic, the minimum inscribed quadrilateral is one which has 
negative infinity for the length of its perimeter. If the inscribed 
quadrilateral be restricted to be non-crossed, the minimum one will 
not be determinable by a general method. 

It will be seen how it is that some sides may have to be consi- 
dered negative, by consideration of figs. 33 and 34, where A' and B' 
are the geometrical images of A and B. In fig. 33, AP + PQ + QB 
( = A'B') is a minimum ; while if we make the same construction in 
fig. 34 we only get AP-f-PQ + QB equal to A'B' and a minimum if 
we consider PQ negative. 

We may also determine the value of the perimeter without con- 
sidering merely a limiting case. 

In fig. 30, let BP = a;, BQ = t/, DR = 2;, DS = w?; PQ = A;, QR = Z, 

RS = m, SP = n; z.BPQ = a, etc. 

Then k = axjosa -f- y cos/?, l = (b - y)cos)8 + (c - z)cosy, 

m — zcosy -H i«?cos8, n = {d- i«?)cos8 + (a - a:)cosa. 

Hence, if p denote the perimeter, 

p^zk + l + m + n 

= acosa + 6cos^ + ccosy -H c?cos8 

= asiny' + 6sin8' + csina' + damft' 

= 2R[sina'siny' + sin^'sinS' + siny'sina + sin5'sin/3'] 

= 4R[sina siny'sin/3'sin6'] 

= (oc -H bdyU = hk/B.. 



26 

General Case. 

In a similar way it may be shown that in any polygon of an odd 
number of sides, there may be inscribed a polygon of the same 
number of sides whose perimeter is a minimum ; and that this in- 
scribed polygon is uniquely determined. The directions of its sides 
and the positions of its vertices may in fact easily be calculated as 
in the case of the triangle. 

In the case of a polygon with an even number of sides no polygon 
of minimum perimeter can be inscribed, unless a relation holds 
among the angles of the given polygon, namely, A-B + C-D+ ... 
must be zero. If this relation holds, an infinite number of polygons 
of minimum perimeter may be inscribed. Any two of these in- 
scribed polygons have equal perimeters and have corresponding sides 
parallel. 

Paeticulab Case. 

The particular case of the regular polygon is of some interest. 

The polygon of minimum perimeter inscribed in a regular odd 
polygon is a regular polygon having its vertices at the middle points 
of the sides of the given polygon. 

Thus, in the case of the pentagon (fig. 35), putting l APT = 

^BPQ = a, etc., we get, successively a = y~€ = ^ ; and therefore 

all the triangles, APT, BQP, etc., are isosceles. Also PB = RD = 

DS = TA = AP ; and thus the vertices of PQRST bisect the sides 

of ABCDE. 

If now we consider a regular hexagon (fig. 36), we do not get 
a = 13, by the reasoning that gave us this result in the case of the 
pentagon ; but only a = y = c ; and ^ = 5 = f . 

But if we suppose a>)8, we have 

BQ>BP, CQ>CR, DS>DR, 
ES>ET, FU>FT, AU>AP; 

whence, adding, we have BC + DE + FA > CD + EF + AB ; which is 
not the case ; and therefore a = f3. 

And in this case P is not necessarily the middle point of AB ; for 
if we make AP = OQ = CR = ES = ET = AU, the hexagon PQRSTU 
will satisfy the conditions of the problem. 

It may easily be shown that the perimeter of PQRSTU is con- 
stant ; it is in fact ^3/2 times that of ABCDE. 
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On some properties of the quadrilateral. 
By R. E. Allardice, M.A. 

§ 1. In a triangle ABC {Gg. 37), BE is made equal to OF ; to 
find the locus of the middle point of EF. 

Take K the middle point of BC and P the middle point of EF, 
then PK is the locus required. For if E' and F' are the middle 
points of BE and CF, the middle point of E'F wiU lie in PK 
(namely, at the middle point of PK) ; and again if BE' and OF' are 
bisected in E" and F", the middle point of E'T" will Ue in PK 
(namely, at the middle point of KB) ; and so on. At any stage we 
may double the parts cut off from BA and CA instead of bisecting 
them. Hence the locus required is such that any part of it, however 
smaU, contains an infinite number of coUinear points j and hence the 
locus is a straight line. 

§ 2. The proof of the above paragraph obviously holds good if 
BE and CF, instead of being taken equal, are taken in a constant 
ratio. 

Hence the proposition may be stated as a property of the quad 
rilateral, as follows : — 

If the points P and Q divide the sides AB and DO of a quad- 
rilateral ABOD in the same (variable) ratio, the locus of the middle 
point of PQ is a straight line. 

§ 3. Generalisation. 

If L and M (fig. 38) divide the sides AB and DC of a quad- 
rilateral ABCD in the same (variable) ratio l:m, then the locus of a 
point P that divides LM in a given ratio is a straight line. 

For this is true, by last paragraph, if we take the ratio LP:LM 
to be 1:2 ; and hence also if we take it to be 1:4 or 1:8 or 3:8 ; or, 
in general, if we take it to be m:2**; and hence it must be true 
generally. 

J^ote. — In this way the whole plane may be divided into a num- 
ber of quadrilaterals whose sides are proportional (but which are not 
similar). 

§ 4. Second Method op Proof. 

Lemma. — Let ABO, DEF (fig. 39) be two straight lines and DA, 
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EB, FC, be perpendicular to AC ; and let further A'B'C (fig. 40) be 
a straight line and D'A', E'B', F'C be perpendiculars to A'C, equal 
respectively to DA, EB, FC ; then if AB:BC = AB'iB'C, the points 
D', E', F', are coUinear ; and conversely. 

The proof of this lemma comes at once on drawing through D 
and D' lines parallel respectively to AC and A'C. 

[This lemma is obviously connected with a property of the 
simlpest kind of homogeneous strain.] 

Now take (fig. 41) BF:FD = CG:GE. Bisect BC and DE in P 
and Q, and join PQ. Draw perpendiculars BH, CH', etc., to PQ. 
Then, obviously, BH = CH', DL = EL'; and HK:KL = H'K':KX' ; 
and hence FK = GK' and FR = RG. 

The more general theorem of (§ 3) may also be proved in this way 

§ 5. Third Method of Proof. 

As neither of the preceding proofs is exactly Euclidian in char- 
acter, it may be as well to add the following proof. 

Let ABC (fig. 42) be any triangle. 

Make AC = AC ; AB' = AB. 

Then BB' is parallel to CC, and P,Q,R,S, the middle points of 
BC, BC, B'C, B'C, are collinear. 

Now BQ = CR = j^ (BA + AC). 

If we make BD = CE, we have still DQ = ER = j^ (DA + AE) ; 
and hence RQ passes through P', the middle point of DE. 

The more general theorem, in which BD:CE is a constant ratio, 
may be proved in much the same way. 

§ 6. By means of the proposition of the preceding paragraphs, a 
simple proof may be given of the following well-known theorem * : — 
The middle points of the diagonals of a complete quadrilateral 
are concurrent. (Fig. 43). 

Make BHDK a parallelogram. 
Then AH/HE = (AH/HB)(HB/HE) 

= (BK/KF)(KC/KB) 
= KC/KF. 
Hence the middle points of HK, AC and EF are collinear. 



* Numerous proofs of this theorem have been giveriy one of the simplest being that 
contained in Taylor^s Ancient and Modern Geometry of Conies, § 107. For some 
account of the history of the theorem, see the last paper in this volume of the Proceed- 
ings, by Dr •/*. S Mackay, 
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§ 7. The following proof, by means of co-ordinates, of the general 
theorem of (§ 3) is so simple, that it may be worth while giving it 
here. 

Put (fig. 38) AL:LB = DM:MC = Xifx ; 

AR:RD = LP;PM = BS:SC =p:q. 

Let the co-ordinates of R, P, Q, be (^i,^i), (^aVi^ (gsj^a) ; ^^^ 
co-ordinates of A be (0^,2/1), etc. 
Then ^^ = {qx^ +px^)/(j) + q), 

^2 = {p(^3 + /^4)/(^ + /^) + qQ^2 + /^)/(A' + /^)}/(« + q) 

= {pQ^ + /^4) + qQ^ + /^i)}/(^ + /*)(p + q\ 
^z=(pxs+qx2)l(p+q)' 

Now we may easily show that if we put 

P = {p{x^ - «3) + q{x^ - x^)}/{X + fx)(p + q\ 

Q = {piVi, - 2/3) + q{yi - 2/2)}/(^ + /*)(p + q\ 

then ^2-^3 = /*P; ^3-^1= -('^ + /x)P; ^i-^, = AP. 

Hence 7/1(^2 - ^3) H- 7/2(^3 - ^1) -H 7/3(^1 - ^2) 

= 7?i/aP - 7y2(A -H /x)P -H ry,XP 

= P{ A(7?, - 7^2) + /a(7/i - -72)} 

= P{A(-,.Q) + /x(AQ)} 

= PQ(-A/a + A/a) = 0. 
Hence R, P and S are collinear. 



An Apparatus of Professor Tait's was exhibited which gives the 
same curve as a glissette, either of a hyperbola or an ellipse. 



Fourth Meeting^ 28^A February 1890. 



R. E. Allardicb, Esq., M. A, Vice-President, in the Chair. 



On the Moduluses of Elasticity of an Elastic Solid according 

to Boscovich's Theory. 
By Sir William Thomson. 

The substance of this paper will be found in the Proceedings of 
the Royal Society of Edinburgh^ Vol. xvi., pp. 693-724 ; and Thom- 
son's Mathematical and Physical Papers^ Vol. iii., Art. xcvii., pp. 
395-498. 
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Fifth Meeting, lith March 1890. 



R. E. Allardicb, Esq., M.A., Vice-President, in the Chair. 



On the different possible non-linear arrangements of eight 

men on a Chess-board. 

By T. B. Spragub, M.A. 

The question having been proposed to me as a puzzle : To arrange 
eight men on a chess-board, so that no two of them shall be in the 
same line, — that is to say, that no two are to be in the same column, 
nor in the same row, nor in the same diagonal line, — I succeeded 
before very long in solving it by finding the annexed arrangement. 
(Fig. 45.) 

Having been subsequently informed that other solutions had been 
found, I set to work systematically to ascertain how many there are, 
with the result that I found that there are twelve essentially dif- 
ferent arrangements, and no more, which satisfy the required condi- 
tion. The number is much larger, if we include the arrangements 
which are really the same, but are presented under a different aspect. 

For instance, the above arrangement may be presented under 
eight aspects (Fig. 47). Here Nos. 2, 3 and 4 are got from No. 1 by 
turning the chess-board round clockwise, so that the side which was 
originally at the top becomes successively the right-hand side, the 
bottom side, and the left-hand side ; each of the four sides of the 
board being uppermost in turn. No. 5 is got from No. 1 by invert- 
ing the position of the men on the board ; and then Nos. 6, 7, 8 are 
got by turning it round as before. 

It is clear that the number of arrangements which satisfy the 
condition that no two men shall be in the same column and no two 
in the same row, is 8.7.6.5.4.3.2; and if we assume that each 
essentially different arrangement (or, for brevity, each distinct 
arrangement) may be presented under eight aspects, we see that 
the total number of distinct arrangements which satisfy the condi- 
tion, is 7.6.5.4.3.2, or 5040. The introduction of the further 
condition, that no two men are to be on the same diagonal line, 
reduces the number of arrangements to 12 ; but this is a conclusion 
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I have arrived at by trial, and I have not been able to prove mathe- 
matically that 12 is the correct number. If we start by putting a 
man on the square marked a in Fig. 46, then, by the conditions 
of the problem, the man in the second column cannot be placed 
on either of the three squares in it marked with a cross ; and the 
man in the third column cannot be placed on either of the three 
squares in it similarly marked. We have thus five squares in the 
second column, on which the man can be placed ; and if it is put on 
the square marked 6, this prevents the man in the third column 
from being placed on either of the three squares marked with a 
circle ; so that there are only two squares in the third column on 
which the man can be placed. If, however, we put the man on the 
square marked c in the second column, we shall find that this leaves 
four squares open in the third column ; while, if we place the man 
on any one of the three remaining squares in the second column, 
there will be three possible positions in the third column. We see 
thus that the number of possible positions in any column depends 
on the squares occupied in the preceding columns, in a manner which 
does not seem to admit of mathematical treatment. 

In order to make an exhaustive examination of all the possible 
arrangements, it is not necessary to examine all the 5040 mentioned 
above. Before describing the process, it will be desirable to adopt 
a notation, by means of which we may indicate any square on the 
board, and any arrangement of men. I number the columns from 
left to right, and the rows from top to bottom, as shown in Fig. 45 ; 
then, in referring to any square, the number of the column is placed 
first and the number of the row, second ; thus the spaces occupied 
by the men in Fig. 45 will be denoted by (1,6), (2,1), (3,5), (4,2), 
(5,8), (6,3), (7,7), (8,4) respectively. When we wish to indicate the 
whole arrangement, it is unnecessary to write down the numbers 
of the columns, as these all run in regular succession; and the 
arrangement in Fig, 45 is sufficiently denoted by (61528374). 
Turning now to Fig. 47, we see that the eight arrangements are 
denoted by the following numbers : — 

(1) 61528374 (5) 47382516 = iXl) 

(2) 57138642 = t;?(l) (6) 42861357 = ir;>(l) 

(3) 52617483 -ir(l) I (7) 38471625 = r(l) 

(4) 75316824 = 77?(1) | (8) 24683175=^^(1) 

On examining these, we observe that (5) is got from (1) by 
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inverting the order of the figures ; and in the same way (6) is got 
from (4), (7) from (3), and (8) from (2). This process I call inver- 
sioii, and I denote it by the initial letter i, so that 

i(61528374) = 47382516. 

If we subtract each of the numbers in (l)from 9, we get 38471625, 
which is (7). This process I call reversion, and denote by r, so that 

r(61 528374) = 38471625. 

Next, inverting the figures, we have 

tr(61528374) = {(38471625) = 52617483, 

which is (3). We have by these processes got the arrangements 
(5), (7), (3), from (1). In order to get the four remaining arrange- 
ments, we must interchange the columns and the rows. Thus in 
(1), (1,6) is changed to (6,1), (2,1) to (1,2), and so on; then arrang- 
ing in the order of the new columns (or the old rows), we get (1,2), 
(2,4), (3,6), (4,8), (5,3), (6,1), (7,7) (8,5), or simply 24683175, which 
is No. (8). This process I call perversion, and denote by p, so that 

jo(61528374) = 24683175. 

Then inversion of (8) gives us (2) ; reversion of (8) gives us (4) ; 
and inversion of (4) gives us (6) ; these processes being symbolically 
denoted as follows : * 

ijt?(6 1528374) = i-(24683175) = 57138642, 

rjt?(61528374) = r(24683175) = 75316824, 

irjt?(61528374) = i(75316824) = 42861357. 

It will now be useful to show how, by means of the numbers, we 

* It is not necessary for our present purposes to investigate the laws 
according to which our symbols of operation, t, r, p, combine with each 
other ; and I will therefore content myself with stating a few of the principal 
laws, without any demonstration : — 

ir= rif ip=pr, rp=pi; 

irp= rip = ipi = rpr =pir = pri. 
As an illustration of the use of these relations, let us take the processes by 
which arrangement (2) in Fig. 47 is got from No. (1) ; that is to say^ th^ pro- 
cess of turning the chessboard clock-wise through a quadrant. We have 
seen that {2)=ip{l), so that the operation will be denoted by ip. If, now, 
we repeat this process, we have {ip)^=ip[ip) = ip.pr=iph=ir; and this, as 
we have seen, is the process by which (3) is got from (1). Again, if we 
repeat the same process once more, we have 

(ip)3 = ip{ip)' = ip.ir = rpr. r = rpr^ = rp ; 
and this, as we have seen, is the process by which (4) was got from (1). 
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may ascertain, with regard to any particular arrangement, without 
any representation of it on a chess-board or in a diagram, whether 
it satisfies the condition that no two men are to be on the same 
diagonal line. For this purpose we need a sort of transformation of 
co-ordinates. If, instead of referring the position of each square to 
two sides of our board, as we have done hitherto, we refer it to one 
diagonal, AB (Fig. 48), and a perpendicular through its extremity. A, 
so as to fix the position of any square by the co-ordinates AM, PM, 
of its middle point, P, taking as our unit the half-diagonal of a 
square ; then a square which is indicated by {xy) according to the 
former plan, will be indicated by (x + y - 1, a; - y), according to the 
latter plan ; or, if we produce the diagonal backwards, and take as 
our new origin a point O, distant half a diagonal from the corner, 
the square formerly denoted by (xy) will be now indicated simply 
by (a + 2/, ic-y). 

According to this notation, the arrangement in Fig. 1 is indi- 
cated by (7,-5), (3,1), (8,-2), (6,2), (13,-3), (9,3), (14,0), (12,4). 
Here the fact that all the eight numbers which stand first in the 
pairs, are different, — 7, 3, 8, 6, 13, 9, 14, 12, — shows that no two 
men are on the same diagonal line perpendicular to AB ; and the 
fact that all the eight numbers which stand second in the pairs, are 
different, -5, 1, — 2, 2, — 3, 3, 0, 4, — shows that no two men stand 
on the same diagonal line parallel to AB. 

We see also that our problem may be stated without any mention 
of a chess-board : Required to arrange two sets of numbers, 1 — 8, 
in eight groups, each containing one out of each set, so that the sums 
of the numbers in the eight groups shall all be different, and the 
differences of those numbers shall also all be different. 



As another illustration, I will show how, by the same operations i, r, p^ 
all the aspects in the text, instead of being got from No. (1), are got from 
one of the others.—say (4). 

We have (4)=rp(l) 

.•.r(4)=r2i>(l)=i)(l), 
andiw(4)=i)«(l) = (l); or(l)=pr(4). 
Then (2) = ipCl ) = ippr{^) = (pM4) = Jr(4) ; 
(3)=ir(l)=irpr(4) = i.ipi(4)=2n(4); 
(6) = i(l) = ii>r(4) = i.ii)(4) =p(4) ; 
(6) = irp(l) = irppr{i) = ir/)V(4) = tr'(4) = i(4) ; 
(7)=r(l)=rpr(4)=n^(4); 
(8)=l?(l)=l^r(4)-r(4). 
It will be noticed that {rp)-^=pr; and similarly (tp)-^=;», («rp)-^=pri=irp. 
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I will now, by means of the preceding notation, explain the pro- 
cess of finding by trial the various possible arrangements. I begin 
by putting a man on the square (1,1) ; then I put a man in the 
second colimm on the highest admissible square, which is (2,3), 
because, by the conditions of the problem, the man must not be 
put on (2,1) or (2,2). In the third column, now, the man must not 
be put on (3,1), (3,2), (3,3), or (3,4) : I therefore put it on (3,5). 
In the fourth column, the square (4,2) is admissible, and I there- 
fore put the man there ; and the man in the fifth column is similarly 
put on (5,4). Looking now at the position of the men on the board, 
we see that it is not possible, according to the conditions of the 
problem, to place a man anywhere in the sixth column. I therefore 
move on the last-placed man to the next admissible square, which is 
(5,8) ; but it is still not possible to place a man in the sixth column. 
I therefore remove the man from the fifth colimm, and move on the 
man in the fourth column to the next admissible square, namely 
(4,7). Then three more men may be placed, thus, (5,2), (6,4), 
(7,6) ; but it is not possible to place a man in the eighth column, 
and the last-placed men have therefore in succession to be moved on, 
and, when necessary, removed from the board. This explanation 
will, I trust, be sufficient to enable my readers, with the chess-board 
before them, to understand the process by which the following 
arrangements were successively arrived at and found inadmissible, 
until the last was arrived at, namely, 15863724. 



13524 X 

8x 

1357246 X 

4x 

1358246 X 

4x 

136275 X 

136824 X 
5 X 

137248 X 
85 X 



13824 X 

7x 

13862 X 

425 X 

14253 X 

8x 

14273 X 
142837 X 

63 X 

1463 X 



1468253 X 

73 X 

146835 X 

7x 

1473625 X 
1473825 X 
14752 X 

82 X 

1483 X 
14852 X 
.3x 



1526374 X 
1528374 X 

63 X 

4x 

157248 X 

63 X 

8x 

1582473 X 

736 X 

15863724 



We can prove that this satisfies the conditions, by transforming it 
as above explained. First writing it in full, we get 1,1 ; 2,5 ; 3,8 ; 4,6 ; 
5,3 ; 6,7 ; 7,2 ; 8,4 ; which for our present purpose may be more 



compactly written 



1 


2 


3 


4 


5 


6 


7 


8 


1 


5 


8 


6 


3 


7 


2 


4 



Writing the transfer- 
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mation iu a similar way, we get 



2 


7 


11 


10 


8 


13 


9 


12 





-3 


-5 


-2 


2 


-1 


5 


4 



and 



since the eight numbers in the upper line are all different, and the 
eight numbers in the lower line are also all different, we see that 
this is an arrangement which satisfies the conditions ; or, briefly, a 
solution. Proceeding thus, we get the following twelve solutions : — 



15863724 . 


. (1) 


25741863 . 


• (5) 


27581463 


. (9) 


16837425 . 


• (2) 


26174835 . 


• (6) 


»35281746 


.(10) 


24683175 . 


• (3) 


26831475 . 


• (7) 


35841726 


.(11) 


25713864 . 


• (4) 


27368514 . 


• (8) 


36258174 


.(12) 



As the work proceeds, various methods of shortening it slightly, 
suggest themselves ; thus, if, when five columns, and consequently 
five rows, are occupied, the remaining three rows are adjacent rows, 
as, for instance, the fiif th, sixth, and seventh ; the problem is reduced 
to placing three men on a small board containing three squares 
in each side ; and it is easy to see that this is impossible under the 
conditions. We therefore know that we cannot complete the 
arrangement, and it is useless to proceed to place a man on the sixth 
column. The same is the case if the two last remaining rows are 
adjacent. 

We have to be on the watch against recording an old solution 
which presents itself under a new aspect. In order to explain this 
more clearly, 1 have given in the following table the numbers indi- 
cating four aspects of each of the twelve solutions. 



(1 

(2 
(3 

(4 
(5 
(6 

!l 

(9 
(10 
(11 
(12 





Reversion, 


Ferversion, 


a 


h — r{a) 


c=p{a) 


15863724 


84136275 


17582463 


16837425 


83162574 


17468253 


24683175 


75316824 


61528374 


25713864 


74286135 


41582736 


25741863 


74258136 


51842736 


26174835 


73825164 


31758246 


26831475 , 


73168524 


51468273 


27368514 


72631485 


71386425 


27581463 


72418536 


51863724 


*35281746 


*64718253 


*53172864 


35841726 


64158273 


57142863 


36258174 


63741825 


63184275 



Reversion of 
Perversion, 
d = rp(a) 
82417536 
82531746 
38471625 
58417263 
48157263 
68241753 
48531726 
28613574 
48136275 
*46827135 
42857136 
36815724 



(1) 
(2) 

(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 

(11) 
(12) 



Here each aspect h is got from the corresponding a by rever- 
sion ; each c is got from a by perversion ; and each d is got from 
c by reversion. It is unnecessary to write down the remaining 
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four aspects of each solution, as they will be got by simple inver- 
sion of the four that are here given ; but it will be convenient to 
denote them by A, B, C, D, respectively, so that A = i(a), etc. This 
order is preferable in some respects to the order adopted in Fig. 
47. We have placed (l)a first because it is the first solution that 
our process gives us. It will be noticed that solution (10) differs 
from all the others, in that we get no new aspects by inverting the 
order of the numbers in it ; thus if we invert the order in (10)a, 
we get 64718253, which is (10)6; similarly the inversion of (10)c 
gives us {10)d, As an example of an old solution presenting itself 
under a new aspect, we may take 17582463, which the above table 
shows us is aspect c of solution (1). Having arrived at this solution 
in the course of our process, we see by inspection of the board that 
it is a new aspect of an old solution ; for such inspection informs us 
that it is an aspect of a solution commencing with 15, which, there- 
fore, we must have got already. It is easy to see further that, when 
in our trials we have got 17 as the commencing numbers, it is use- 
less to make trial of arrangements that have a 2 in either the 3rd, 
4th, 5th, or 6th place ; because, if any one of these should give us a 
solution, it must be another aspect of one that we have already got. 
For a similar reason, if we begin with 14, it is useless to make 
trial of 142 . . . . ; if we begin with 15, it is useless to make trial 
of 152 . . . . or of any arrangement containing 2 in the 4th place ; 
if we begin with 16, we similarly reject arrangements with 2 
in the 3rd, 4th, or 5th place ; and when we begin with 18, we reject 
arrangements which have 2 in the 3rd, 4th, 5th, 6th, or 7th place. 

There is no similar proposition when we begin with 2 ; but when 
we begin with 3, it is useless to try arrangements beginning with 
31 ; for by interchanging the columns and the rows, each of those 
arrangements must give us an arrangement commencing with 2 ; 
and all such we have 'already had. Similarly, if we began with 4, it 
would be useless to try any arrangement which has 1 in the 2nd or 
3rd column ; if we began with 5, it would ba useless to try any 
arrangement which has 1 in the 2nd, 3rd, or 4th column ; and so on. 

The new aspects which we have thus far considered, are those 
which are contained in the column c, and are obtained by the process 
I have termed perversion. In order to recognize the others, we must 
first see whether there is a man nearer to a corner square than the 
one we started with. For instance, our process leads us to the solu- 
tion 35714286 ; and as the square (7,8) is next to a corner, this is 
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a solution we have already got, and our table shows us that it is the 
D aspect of No. (6). Even if a solution contains no man nearer to 
a comer than the one we started with, it may still be a new aspect 
of an old solution : for instance, 36824175 contains 2 other men 
that are at the same distance from a comer as the one we started 
with, namely, in the third square ; and the arrangements commenc- 
ing with these respectively are 35841726, which is No. (11); and 
37285146, which is the B aspect of the same; while 36824175 is its 
C aspect. When we have completed all the trials that begin with 3, 
it is unnecessary to proceed any further ; for, putting a man on the ' 
square (1,4), we see at once that, if no man is to be nearer to a comer 
than this man is, the only admissible arrangement of the four men 
on the extreme columns and rows is the one shown in Fig. 49 ; and 
on trial we find that only one other man can be placed on the board 
consistently with the conditions of the problem. 

As regards the remaining four squares in the first column, we 
see that by commencing with any one of them we shall simply get 
the h aspects of the solutions we have already got; for instance, 
if we put a man on the 7 th square, we shall get the h aspects 
of the solutions we got by putting the first man on the 2nd 
square. 

On examining the twelve solutions we have got, the first thing 
that strikes us is the general absence of anything like symmetry, or 
apparent law, in them. It is, of course, clear that the conditions of 
the problem preclude bilateral symmetry ; but we may have centric 
symmetry, so that to every man on the board there corresponds 
another at the same distance from the centre on the other side of it. 
One of our solutions is of this character, namely No. (10). In the 
numbers indicating it, 35281746, and 53172864, we observe that the 
sum of each pair of digits equidistant from the two ends, is 9 ; the 1st 
and the 8th, the 2nd and the 7th, and so on. It results from this 
centric symmetry that, as already mentioned, there are only 4 aspects 
of this solution, and not 8, as there are of the others. Hence the 
total number of arrangements which satisfy the conditions of the 
problem is 8x11+4 = 92; the total number of possible arrange- 
ments that contain one in each column and one in each row, being 
8.7.6.5.4.3.2 = 40,320. 

From a study, either of the board or of the numbers in our table, 
we see that some of our solutions can be deduced the one from the 
other ; for instance, if in No. (3)a we remove 2 from the first place 
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and put it in the last, we get 46831752,* which represents No. 
(4)A. In order to see whether a transformation of this kind is 
possible we may proceed as follows, referring the men to a diagonal 
of the board and its perpendicular, as pre\'iously explained. 



Columns x 


.-1, 





Rows y 


. 7, 


5 


x + y 


. 6, 


5 


x-y 


.-8,- 


-5 



1, 2, 3, 4,5,6, 7, 8 

2, 4, 6, 8, 3, 1, 7, 5 

3, 6, 9, 12, 8, 7, 14, 13 
-1,-2,-3,-4,2,5, 0, 3 



9, 10 

2, 4 

11, 14 

7, 6 



When we remove 2 from the first place to the last, and thus get 
46831752, the new values of x + y smd x-y are 11 and 7, which 
are difierent from any of the others ; and this shows that the new 
arrangement is a solution. So, if we take 5 from the last place to 
the first,! ^^ S^^ * solution 52468317, which is (8)0. But if we 
attempt to repeat the process in either direction, we get no new 
solution. Taking, for instance, 75246831, we have a value oi x + y, 
namely 6, which is the same as one of those we have already got ; 
and this shows that the arrangement is not a solution. Similarly, if 
we take the arrangement 68317524, we get a value oi x + y, namely 
14, which we have already got ; and this arrangement also is not a 
solution. I indicate this connection between the three solutions 
graphically as follows : — 

X (8)0 (3)a (4)A x 

Again, if we take the solution (ll)c, namely 57142863, and sub- 
tract 1 from each of the numbers except 1, which is to be replaced 
by 8, we get 46831752, which is (4) A. This operation I denote by 
8 (the initial letter of subtract), so that 46831752 = s(57142863), and 
(4a) = w(llc). Each of these operations, s and tj has a simple inter- 
pretation with reference to the arrangement of the men on the chess- 
board, which is so easily understood when the men are placed on the 
board, that it is unnecessary to explain it. I represent (l^)c 

the connection between (ll)c and (4) A graphically thus : — 

(4)A 

It will be found that seven of our solutions, namely, Nos. (1), 

(2), (3), (4), (6), (8), and (11), can be deduced one from the other by 

r 

* This operation I denote by t (the initial letter of transpose), bo that 
46831762 = t(24683176); whence (4A) = «(3a). 

t This operation will be denoted by t-\ so that 62468317 =t-i(24688 176) 
and we have (4A)=t(3a)=t3(8C). 
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combinations of the s and t operations ; and I think the reader with 
the board before him will have no difficulty in understanding, with- 
out further explanation, the following diagram, which represents the 
connection between them. 



(8)0— (3)a— (4)A— x 



X 



X 



X 



X X 



X 



X 



X 



X 



X 



X 



X 



X 



X- 



X 



X 



X — (2)d X 



X _(8)D — (3)5— (4)6 X 



X _(1)C— (6)0— (ll)c^— X 



X 



(1)D— (6)D— (ll)c— X 



(8)0— (3)a— (4)A x 



(8)0 



X — (2)c X 



X - (1)D 



.(8)0 



It will easily be seen that 

(80) = s(lD) = 8H(2c) = sH\Ud) = sH%ib) 

= s¥(6C) 
= etc. 

The remaining 5 solutions, Nos. (5), (7), (9), (10), and (12), do 
not admit of being transformed in a similar way. It will be noticed 
that (2)c and (2)d are derived from the adjacent arrangements in a 
different way from the rest, namely, by the compound operation st, 
or 8~H~^, or 8t~\ or 8~H. Thus (2)c is got from (1)D by carrying a 
man from the right hand bottom corner of the board to the left hand 
top corner, and moving all the other men diagonally through one 
square downwards to the right. These two solutions have each a 
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man in a comer square of the board, and I propose to call them 
" comer solutions ". 

Although I have satisfied myself by trial that there are 12 solu- 
tions and no more, I have not been able to discover any means of 
proving otherwise that this is the correct number. In order, if pos- 
sible, to discover a law for the number of solutions, I have investi- 
gated the problem for the cases where the board has 4, 5, 6, 7, and 9 
squares in its side ; and the results are shown in the following table. 



Number of 


Distinct 


Of which 


Total 


Sides. 


Solutions. 


are Centric. 


Solutions. 


(1) 


(2) 


(3) 


(4) 


4 


1 


1 


2 


5 


2 


1 


10 


6 


1 


1 


4 


7 


6 


2 


40 


8 


12 


1 


92 


9 


46 


4 


352 



Total 
Arrangements. 

(5) 



24 

120 

720 

5,040 

40,320 

362,880 



Quotient 
(5)--(4) 




There is no law visible in the progression of these numbers. 

The solutions in the case of the boards containing 4, 5, 6, and 7 
squares in a side, are given in Fig. 50. 

The solutions (1), (2), (4), (6) of the 7-board can be derived one 
from the other in an endless series, as shown in the following scheme. 
Here I have not thought it necessary to indicate the different aspects 
of each solution, but the 1 printed in heavy type is a different as- 
pect from the 1 printed in ordinary type. 

1_6— 6— 1—2—4— 2— 1—6— 6— 1-2-4— 2—1 



1—2—4—2—1—6—6—1—2—4—2—1—6—6-1 

I I I I I I I I I I I I I I I 
2—1—6—6-1— 2—4—2—1—6—6 _1— 2— 4— 2 

The succession of the numbers vertically is easily seen to \ 
1_2— 6— 4— 6— 2— 1— 1. 

The other solutions (3) and (5) cannot be similarly derived frc 



- nnv other. 



c^_j.i,„ 1 1 u:^u 1 a 
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136824975 
—7285946 
—8692574 
146392857 
— 825397 
—7382596 
— 925863 
—8397526 
157938246 

42863 

—9642837 
168374295 
174835926 
^9625 



II. (1) 

(2) 
(3) 
(4) 
(5) 
(6) 
(7) 
(8) 
(9) 
(10) 

(11) 
(12) 

(13) 

(14) 

(15) 

(16) 

(17) 
(18) 
(19) 
(20) 
(21) 
(22) 



241796358 
—7139685 
—8396157 
—9731685 

53168* 

257936418 

48136 

—8136974 

96374 

— 693147 

^74 

—9418637 

261379485 

— 753948* 

— 958473 

—3184975 

-9358417 

275194683 

—9631485 

281479635 

—5396417 

—6931475 



III. (1) 358296174 
^7146 



I 



3) —9247186 
[4) 362951847* 
;5) —8159247* 
6) — 519724 

(7) —9741825 

(8) 372859164 

(9) 386192574 



IV. (1) 427918536 



Of these solutions, the four marked with a star are centric, 
namely 11.(5), 11.(14), 111.(4), and 111.(5) ; and it may be useful 
to give the graphic representations of them (v. Fig. 51). 

In our search for a law, an obvious idea is to try whether any 
relation can be found to exist between the numbers of the solutions 
when the board has n squares in its side and when it has n+1. 
When we examine from this point of view the solutions we have 
found, we immediately see that it is only in exceptional cases that 
a solution in the one case can be deduced from one in the other case. 
If we take a solution for the board containing (n+ 1) squares in its 
side, which for brevity we may call an (n + l)-board, and remove from 
the board one of the outside columns and one of the outside rows, 
we shall get an n-board, but we shall not in general get a solution ; 
for, as there is a man on each column, and a man on each row, we 
shall have removed two men, unless it happens that there is a man 
on the corner square of the board which is common to the column 
and the row which we have removed. Conversely, when a solution 
for the n-board contains no man on one of the diagonals, we can get 
from this solution two solutions for the (w+ l)-board, by producing 
the free diagonal, and adding a new column and a new row so as to 
intersect in the prolongation of the diagonal, first at the one end, 
and then at the other, and placing a new man on the corner square 
thus obtained. In this way solution (1) for the 5-board is got from 
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3 solution for the 4-board ; solutions (1) and (2) for the 7-board are 
<t from the solution for the 6-board ; solutions (1) and (2) for the 
board are got from solution (3) for the 7-board ; and the following 
3lutions for the 9-board are got from solutions for the 8-board : — 

Nos. 1, 6, 10, 14 for the 9-board from No. 4 for the 8-board. 

7, 13, „ „ 5 

2, 5, „ „ 6 

8, 11, ,, „ 7 

3, 12, „ „ 9 

4, 9, „ „ 10 

By examining the diagrams for the solutions for the 9-board, it 
is easy to see that they will give us 32 solutions for the 10-board. 
The numbers of these " corner solutions " are — 

For the 4-board, none 

>) 
)) 

Here, again, there is no obvious law in the series of numbers. 

It is sometimes possible to get a solution for the n-board from 
one for the (w + 1) board, by removing a man from the board, and 
both the column and the row containing it, and then closing up the 
board ; for instance, solution (2) for the 5-board gives in this way 
the solution for the 4-board. Conversely, one or more solutions for 
the (n + l)-board may sometimes be got by inserting a new column 
and a new row, and placing a man on the square in which they inter- 
sect. In this way we get : — 

From Solution (1) for the 8-square, 1 solution for the 9-square. 

(2) .. 1 



5 




1 


6 




none 


7 




2 


8 




2 


9 






14 
32 



)) 



,) 



>J 



)> 



JJ 



>> 



J> 



if 



>» 



)> 



>> 



(3) 





(4) 


5 


(5) 


3 


(6) 


2 


(7) 


4 


(8) 





(9) 


3 


(10) 


3 


(11) 


1 


(12) 


1 


Tot 


AL, 24 



)i 



>> 



>> 



>J 



>> 



)J 



n 



>» 



>> 



)) 



»» 
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In consequence, however, of two of these solutions being given 
in 2 ways each, the actual solutions got are only 22, or 8 in addition 
to the comer solutions. There are thus 24 solutions for the 9-board 
which cannot be got from solutions for the 8-board. This method 
of treating the question therefore does not lead us to a law. 

A study of the diagrams suggests several interesting propositions 
for investigation ; for instance, in no one of the solutions we have 
got for the 6-, 7-, 8-, and 9-boards, is there an arrangement of 4 men 
such as we have them in the 4-board solution ; and I am inclined to 
think it is impossible there should be, but I have not succeeded in 
proving it. Similarly, I am inclined to think it is impossible there 
should be an arrangement of 4 men on the outside columns and 
rows of a board, as shown in Fig. 52 ; but this also I have not 
succeeded in proving. 



On the equations of Vortex motion, with special reference 

to the use of polar co-ordinates. 

By C. Chiiee, M.A., King's College, Cambridge. 

In several previous communications* to the Society, I have con- 
sidered the equations of vortex motion in two dimensions in a com- 
pressible fluid. In the present communication I propose to consider 
certain forms of the hydro-dynamical equations of a more general 
kind. In certain cases the fluid will be supposed to be rotating, 
prior to the introduction of the vortex motion, with uniform angular 
velocity about a fixed axis. 

Using the same notation as in my previous papers, and supposing 
the axes of x and y rotating with uniform angular velocity w about 
the axis of z, which is supposed fixed, we can easily prove by the 
method of my paper " On vortex motion in a rotating Fluid " f that 
the equations of three dimensional motion are the following — 

^w O 2 ^ ^P ^^ /1\ 

ot p ax 

^ + 2<oM-<oV=--J + Y (2), 

0^ p ay 

^^=-JLf + Z (3). 

Sc p dz ^ ^ 



Proceedings, vol. V., p. 52 ; vol. VI., p. 59 ; vol. VII., p. 29. 
t Proceedings, vol. VII., p. 29. 
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The component velocities u and v are velocities relative to the 
moving axes. It will also be remembered that 

8 d d d d / iv 

denotes differentiation following the fluid ; while — denotes variation 

dt 

at a point fixed relative to the moving axes. 

The equation of continuity is simply 

1 8p du dv dw _.. /^. 

7¥ + ^ + 5y + ^-" *^^- 

The ybrm of this equation is, it will be noticed, the same whether 
(0 exist or not. 

Employing the usual notation for the components of the appa/rent* 
vorticity — 

f = ^l — - _ — I, etc., ... ... ... (by, 

we easily deduce from the preceding equations in the case when the 
components X, Y, Z of the external forces vanish or are derivable 
from a potential, the following equations f for the variation of the 
vorticity — 

.Km \^'^ +.,^ ^,y^ x^^l m 



'h 



S(r,/p) 



P-^ 



8(f +<-)//» 



St 






8(j) 



In the left-hand side of (9) we have introduced which is of 

8t 

course zero. The components ^, rj are of course absolute as well as 

apparent vorticity components along the instantaneous directions 

of the corresponding axes. These equations will remain unaltered 

when we transfer the origin to any point fixed relative to the original 

rotating axes, the new axes of x and y being supposed to rotate with 

the same angular velocity w as the original. In proof of this it is 

sufficient to remark that the co-ordinates of any point, whether fixed 

or moving, referred to the new axes will differ from those referred to 



* See Proceedings, vol. VII., p. 32. 
t Lamb's Motion of Fluids, Note D. 
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the old only by constant quantities independent of the time, and so 
all velocity and vorticity components are unaffected by the change. 
The components of velocity and vorticity must satisfy (7), (8), 

(9) and in addition the ordinary boundary conditions and the equa- 
tion), (5), of continuity. This last equation requires inter alia that 
two contiguous elements of fluid the one inside and the other outside 
a vortex filament have the same component of velocity along the 
normal to the surface of the filament. 

As the fluid is supposed frictionless it is not necessary, unless it 
be involved in the equations (7) - (9), for the component velocities 
in elements just inside and just outside the surface of a vortex fila- 
ment to be identical in any direction which lies in the tangent plane 
to the surface of the filament. 

Suppose now that when cu is zero the components of velocity and 
vorticity in the fluid constituting a vortex ring, whose centre is in 
the axis of z and whose plane is parallel to xy, are given by 

u = vcosO v = vsinO w = Wq\ .^x^^ 

^=-Qsin^ 7/ = Qcos^ f=0 / '•• ^^"/' 

where = \^n~\yjx)^ and v, Wq and 12 are independent of 9, 

These components will thus satisfy (7), (8), (9) and other neces- 
sary conditions when cd is zero. When, however, cd ceases to be zero 
the equations (7), (8), and (9) are no longer satisfied by (10). These 
equations may, however, be all satisfied by adding to the components 

(10) the additional terms 

u = (orsin^ v= - (orcosO t^; =» \ .^ ^ v 

^=0 7/=o f=-W ••* ^ ^' 

where 7^ = o(P + y\ 

The additional terms in the velocity components add nothing to 
the components ^ and rj of vorticity, and they are consistent with 

(= -0). 

To prove our statement we must substitute in (7), (8), and (9) 
the system of velocities and vorticities obtained by combining (10) 
and (11). Doing so and remembering that (9) is satisfied by the 
values (10) when (o = 0, we find that it is still identically satisfied, 
and that (7) and (8) respectively lead to 

- sin^ — - _ ilcosO — - + — ?^?_ -^ = — .--(vcos^ + wrsin^) ... (12), 
ot ot p ht T du 

cos^ — - - 12sin^ -— i: = — -— (vsin^ - wrcos^) ... (1 3). 

0^ ot p ot T du 
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Multiplying (12) by cos^, and (13) by sin^, then adding and 
dividing out by 12 we simply get 

SO 



= — Ct). 



8e 

This signifies that the fluid in the vortex ring moves relative to 
the rotating axes with angular velocity - w about the axis of the 
ring, which is exactly the motion indicated by (11). 

Again subtracting (12) multiplied by sin^ from (13) multiplied 
by cos^, we get 

Sn_ il 8p ^nv 
St p 8t r 

which is precisely the same relation as when w is zero. 

The additional terms (11) in the velocity of the fluid possessed of 
vorticity have their resultant in the tangent plane at every point on 
the surface of the vortex ring, and so their existence in no way 
affects the equality of the velocity components inside and outside 
the ring in the direction of the normal to the surface. 

This indicates that a circular vortex ring may exist in a rotating 
fluid with its plane perpendicular to the axis of rotation, provided it 
have, in addition to the motion existing in a similar ring in a non- 
rotating fluid, a uniform angular velocity about the axis of the ring 
which is equal in magnitude but opposite in direction to that of the 
undisturbed rotating fluid. As in the case of a non-rotating fluid, 
there may be an infinite plane boundary at right angles to the axis 
of rotation, the necessary conditions being satisfied by the existence 
of an image on the remote side of the plane. 

If, for instance, we imagine the state of matters in the earth's 
atmosphere in latitude A the same as in fluid limited by an infinite 
plane coincident with the tangent plane to the earth's surface, the 
fluid rotating with angular velocity (osinA about an axis answering 
to the vertical at the point considered ; then a vortex ring whose 
plane was horizontal would have an angular motion wsinA in azimuth 
from east to west through south. In other words, a diameter of the 
ring connecting two definite material cross sections would appear to 
an observer at the point to rotate about the vertical in precisely the 
same way and at the same rate as the plane of vibration of a Fou- 
cault's pendulum. 

Hitherto I have said nothing as to the motion in the fluid sur- 
rounding the vortex ring. When we consider this point a difficulty 



the ring, Thus -p-, -rr- and -;- do not vanish in the fluid surround- 
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appears which I do not see my way to answer satifactorilj. It 
seems well worthy of notice. 

The equations (7), (8), and (9) apply at every point in the fluid, 
and so must be satisfied by the fluid surrounding the vortex ring we 
have just considered. Now the velocity components due to the 
action of a vortex ring in the surrounding fluid satisfy the conditions 
of irrotational motion, but they depend both in magnitude and direc- 
tion on the distance from the instantaneous position of the plane of 

du dv J dw 
— -, — and — 

dz dz dz 

ing the ring, and so in accordance with (7), (8) and (9) unless f + w 
vanish everywhere, the whole of the fluid surrounding a vortex 
ring will take up some form of vorticity of a complicated char- 
acter. 

Theoretically as — , etc., do not absolutely vanish, however dis- 

dz 

tant the point considered from the vortex ring may be, in order to 

avoid setting up vorticity possessed of components ^, ?/ we require 

some species of motion set up whereby {— - w all through the fluid. 

This would require either the motion of the rotating fluid entirely 

to stop, or in addition to the rotation cd about the original axis a 

rotation - w about some parallel axis. The two rotations would give 

a uniform translational velocity the same at every point occupied by 

the fluid. 

In practice when we consider that every fluid is possessed of 

viscosity, and that the terms — etc., due to the action of an ordinary 

dz 

vortex ring become insensible at a small distance from the vortex, 
the conclusion we seem led to is that the vortex ring itself and the 
immediately adjacent fluid will possess, relative to the moving axes, 
a rotation about the axis of the ring in the direction opposite to that 
in which the undisturbed fluid rotates. Throughout the ring itself, 
if its vorticity be great, this angular velocity may approach indefi- 
nitely near - w, but it will gradually diminish in the external fluid 
as the distance from the ring increases and altogether vanish at dis- 
tances where the direct action of the ring becomes insensible. 

The previous remarks also indicate that any variation in any 
velocity component with the distance from a fixed plane perpendicu- 
lar to the axis of rotation, whatever may be the exciting cause. 
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necessitates the existence of some form or other of apparent vortex 
motion. 

For certain purposes it is desirable to employ polar co-ordinates 
in the hydrodynamical equations. Such equations are obtained by 
Basset in his Treatise on Hydrodynamics* His proof, however, is 
not of a very elementary character, and I think the following method 
shows more clearly the meaning to be attached to the symbols 
employed. For some of the results I shall refer the reader to a 
previous paper in the Proceediiigs.\ 

As in the paper referred to, the element of volume has three of 
its edges 8r, rSB and rsin^S<^ intersecting in P, the point whose co- 
ordinates referred to a fixed point 0, and the ordinary polar system 
are r, ^, <^. The first of those three line elements and the tangents 
to the other two at P form a system of orthogonal axes called the 
fundamental axes at P. If, as in the paper referred to, the comer 
of the element of volume opposite to P be termed S', the co-ordinates 
of S' are r + Sr, + 80, + 8ct>. 

Also neglecting squares of the small quantities 80 , 8<^, the cosines 

of the angles between the fundamental axes at S' and those at P 

are given by the following scheme, reproduced from my previous 

paper, — 

atS' 



I 



A 

r <f> 


1 


-80 


- sin08(t> 


80 


1 


- cos08<l> 


siD.08(fi 


cos^8<^ 


1 



r 

titV < 

Let us first suppose that the axes are fixed, so as to avoid dealing 
with too many difficulties at once. 

Let the velocity of a given element of fluid at the time t at the 
point P have the components u, v, w along the fundamental axes at 
P. After a short interval t let this element be at S', so that 

8r = UTj 80 = T, S<^ = — r— ^T. 

r rsinO 

Then the velocity components at time t + r along the fundamental 

axes at S' are 

, du . dues . du^n , dUes, . 

Ui = u+ — T -I- — or + —pfiO + — o</)., etc. 
dt dr dO dcfi 



* Vol. II., Art. 470. 



t Vol. III., p. 109. 



49 

Thus, referring to the above scheme, the components of the velocity 
of the element at time t-\-T, along the fundamental axes at P are 

Wi — v^W - WjSm68(f), etc. 

But the change per unit time in the component of the velocity in 
any given direction equals the accelerating force per unit mass in 
that direction. If then U, V, W be the components of the external 
forces per unit mass at the point P, proceeding to the limit when t 
vanishes, we easily deduce for the equations of motion — 

^«-':'t«''=-l^+u (14), 

0^ r p dr 

T7 + — -— cot<9= -— -^ + V (15), 

ot r r prdu 

8w , icw , vw .^ \ dp , ,,7. /,/jv 

— + _+— cotl9=~_^^^ ^ + W ... (16); 

ot r r fyrsmu d<l> 

where S d , d v d w d ,.„. 

—- = —- + u—- ■\- - + — ... (17) 

6< dt dr r dO rsinedcfy ^ ^ 

signifies as usual differentiation following the fluid. 

These equations agree with Basset's, except that he has u instead 
of w in the last term on the left-hand side of (16). I believe, how- 
ever, that the above is the correct equation. 

By employing the method of flux it is very easily proved that the 
equation of continuity is 

1 8p 1 diwi"^) 1 C? / . /jv , 1 dw ri /^Q\ 

p ot IT dr rsin^ dO r^inO d<\> 

Let us next suppose that the entire system of fluid and axes is 
rotating with uniform angular velocity w about the fixed axis ^ = 0. 

Let u\ v\ w' denote the component velocities at time t at the 
point P (r, ^, </)) relative to the fundamental axes there, which are 
supposed fixed relative to the rotating axes ; and let u, v, w be the 
absolute velocities in the direction of lines fixed in space with which 
the fundamental axes at P coincide at the time t. Then 

u = u\ v = v\ w — V) -{■ {iirsind ... ... (19). 

Let R represent the position at time ^ + t, where t is very small, 
of a point rigidly connected with the moving axes which at the time 
t coincided with the point P ; then the components of the velocity 
at R at time t + r relative to the fundamental axes there, which 
move with the fluid, are respectively 
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, . du , , dv' , . dw 
u + T, V + — T, w + T ; 

dt dt dt 

and so the component velocities at R at time t + T relative to fixed 
axes, coinciding with the instantaneous position of the fundamental 

axes there, are 

, , du' , , dv' . dw' , . /J fCKw 

u -^ — T, V -\- — T, w -{■ T + wrsmC^ . . . (20). 

dt dt dt 

For the last result it must be noticed that referred to axes fixed 
in space, coinciding with the position of the axes at at time t, the 
r and 6 co-ordinates of R and P are the same. From this and the 
consideration that the <^ co-ordinate of R relative to the fixed axes 
exceeds that of P by (ot^ it also follows that the components (20) 
are identical with 

, du , du . dv , dv dw , dw /oi\ 

W + T + (DT V + T + (J)T . W+ T-hiOT ... (21). 

dt d<f>' dt d<t> dt d<t> 

Thus, comparing (20) and (21) and remembering (19), we find 

du' du . du dv' dv . dv dw' dw . dw 
dt dt d(t>' dt dt d<t>' dt dt d<f> 

But from (19) 

du du' dv dv' dw dw' 

d<f} d(f> ' d(f> d<f> d<f) d<f>' 
and therefore 

du _ du du' dv _ dv dv' dw _ dw' dw' 

dt dt diji dt dt dcji dt dt dij> 

Substituting these values and the values 

du du' dv dv' dw dw' 



dr dr^ dr dr dr dr 

du dw dv dv' dw dw' 



+ cusin^. 



+ (orcos^, 



dd dO' dO dO' dO dO 

in the equations (14), (15) and (16) we get after reduction, dropping 
the dashes so that in the following equations w, v, w are the velocities 
relative to the tnoving axes, 

Su v\ (w + o>rsin6f __ 1 rfp^jT (22), 

8t r T p dr 

8v , uv / , . /jv2Cot6/ 1 dp ^y^ /oQ\ 

— + — - (w + iorsmdr =- ^^ + v ... ... (23), 

8t r ^ ' r prdO ^ ^ 

^ + ^ + ^!!^^%2W^.sin^ + i;cos^)=-l-L^^ + W...(24). 
8t r r p r&iii6d(t> 
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In these equations -— has the same form as in (17), and it also 

ot 

as there signifies differentiation following the fluid. The reader will, 

I think, find no great difficulty in this proof if he clearly realise that 

-—— in (20) signifies the rate of change with the time of the com- 
dt 

ponent of the velocity relative to a certain moving axis, at a point 

whose co-ordinates are fixed relative to the moving axes ; while 

-— in (21) signifies the rate of change of the component of the velocity, 
dt 

relative to a certain axis fixed in space, at a point which is absolutely 

fixed in space. 

As a simple illustration of these formulae, let us apply them to 
the case of the earth's atmosphere regarded as surrounding a rotating 
sphere of radius R, the acceleration due to whose attraction at dis- 
tance r from its centre is gW/r^, 

For the undisturbed condition, when the atmosphere is at rest 
relative to the sphere, we have u = v = io = 0. Substituting these 
values along with U = - g^^/r^i V = W = 0, we reduce the equations 
(22) - (24) to 

^ = (oVsin^^ — g^'jir 

p dr 

— -^ = a)Vsin^cos^ V (25); 

p du 

1^ = 

p dcfi 

whose solution is, according siS p = kp or p = kpy, 



log(p//>o) = ^^sin=^-f-R(l-^) 



w / p\^ r (26). 

°' py-^ - p„T-> = Lzi{a,Vsin=e - 2<,r(i - ^)} 

In both results pQ represents the density of the atmosphere at the 
surface of the earth at the poles. 

These results may of coui'se be obtained in a more elementary 
manner by the introduction of the idea of " centrifugal force," as in 
Besant's " Hydromechanics" 

In a paper* already referred to, I obtained the following expres- 

* Proceedings^ Vol. III., p. 113. Cf. Basset's Treatise on Hydrodynamics^ 
Vol. I., Art. 18. 
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sions for the components of vorticity at the point (r, ^, </>) with 
respect to the f undamentiil axes there, 

^ = _1 /_L^-l(^.r)| V ... (27). 

^ 2r \dr ' do) 

When Uj v, w are velocities referred to iixed axes, the above are 
components of the absolute vorticity, but when the velocities are 
relative to moving axes they are components of what I have called 
the ajyjyarent vorticity. 

In order to apply our polar equations to vortex motion we 
require to obtain from (22), (23), (24) a series of equations corre- 
sponding to (7), (8), (9), when as there the external forces vanish or 
are derivable from a potential. Tlie necessary algebraic operations 
are somewhat tedious and lengthy, so I shall merely indicate the 
method of procedure and give the results. The verificatuyii of the 
results will present no great difficulty to any one who keeps their 
form in view in grouping his terms, and who does not allow the 
length of his intermediate expressions to alarm him. 

Replace in every case -— by its equivalent in (17) and notice that, 

ot 

as the external forces have a potential, 

dr r do rsinO dff) ' 

where F is some function of r, 0, <^. 

To obtain the tirst equation multiply (24) by sinO, then differen- 
tiate with respect to 0, and from the result subtract (23) after 
differentiating it with respect to </>. Finally multiply up by r/2sm0. 

To obtain the second equation multiply (22) by cosec^, then 
differentiate with respect to <^, and from the result subtract the 
result obtained by multiplying (24) by r and then differentiating it 
with respect to r. Finally divide out by 2?^. 

To obtain the third e<iuation multiply (23) by ?•, then differentiate 
with respect to ?', and from the result subtract (22) differentiated 
with respect to 0. Finally divide out by 2r^sin^. 

In each case one of the objects of the preliminary operations is 
to eliminate ;; and the external forces. 
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In putting the resulting expressions into the following concise 
forms the identity (17), the equation of continuity (18), and the 
identities (27) are alone required. We thus finally obtain 

ot\ p f L dr r du rsinO a4>j 

- ^%r\ —^rainOlnsinO + rcos^) 4- rcosf/— \ ... (28), 
sin^ \dO d<fiJ 

8t\ rp f L dr r dO r^mdd(fi\ \rj 

4- -^-—-[ —rsinOiusinO + 2;cos6') + sin6^ '-^'l . . .(29), 

8t\r^V[Spf L dr r dO rsin6^G^</)J Vrsint// 

r n d sin^ ^1 / ^" \ /oA\ 

+ CU cos(9_- l_^_^) (30). 

L dr r duj Xrsinc// 

So long as w is not zero these equations apply to the apparent 
not the absolute vorticities. 

In the ordinary case when w is zero and the axes are fixed, the 
following very concise form of the equations seems worth recording. 
Let 12 denote the resultant vorticity, now absolute, and ds an element 
in the direction of the axis of the resultant vorticity at any point in 
the fluid, then it is obvious that 

^± + ±£+J—£^il± (31); 

dr r dO i'ninO d(f> ds 

and so the equations of vortex motion may be written 

plLJ\=nA(^ ... (30a). 

^ St\rsmOpf dsXrsinOf ^ ^ 

In employing polar co-ordinates it will also be convenient to 
possess equations answering to those in Cartesian co-ordinates of 
which in Lamb's* notation the type is 

(fP d^ dM 

u= — + — 

dx dy dz 



Treatise on the Motion of Fluids^ Art. 129. 
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To transfer to polar co-ordinates let as usual 

ar' + y* + 2;^ = 7^, zjr = cos^, yjx = tan<j!>. 

Then dr - a j^ dS cos^os</> d<i> __ sin</> . 

— - — sint7COS9, — , — — •- — ; — Tj., etc. 

dx dx r dx rmnu 

Using these and taking the resultant along the three fundamental 
axes at each point of the components of the velocity in Cartesians, 
and denoting now the velocities along the fundamental axes by 
Uj V, Wy we easily find 

u = — -f — I cos</) -— - sind)cot0 — I M 
dr rl ^dO ^ d<l>i 

ir..rf, , .^ dlj \ dN 
rl ^dO ^ rf</>J r f/</>* 

1 c?P . r . ^ d . cos</) c^ 1 T r j^ d sin<f> rf 1 ■•*■ 
V = . + 1 sm0 — + — — -_ I L - I cos<^ — - — T I M 

r du L dr r a</)J L dr r d<f)j 

cot6 dN 
r dfjy 

tv^J^'^ + Fcos^ ^ _ sin^ rf 1 fLe^g^ ^^^^^1 
rsmed4> L dr r dOM ^ ^J 

[. /, rf . COS^ C? 1 XT 
dr r dO] 

Putting 

R = sin^(Lcos<^ + Msin</)) 4- Ncos^ '\ 

S = cos(9(Lcos</> + Msinc^) - Nsinl9 V ... (32 

T = Mcos</) - Lsinc^ 3 

the above equations reduce to 

u = ^+ ^ /4(Trsin0) - -(Sr)| 



d<t> 

r dO rsin^\ (/</> 

1 dF 1 (d{Sr)_dR\ 
rsinO defy r \ dr dO J 



I dP ^ 1 (dR C^/m .„m\ 



W 



If for shortness we represent the element of volume by dY 
employ the equation of continuity we can replace equation (1 
LamVs Art. 129 by 

P^lffflVl^. 

47rJ ]} p' 8t r, 

a form applying to polar as well as Cartesian co-ordinates. 
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The integration is to be extended to every portion of space 
wherein the fluid is varying in density, and r^ denotes the distance 
of the element dV' wherein the density is p' from the point where 
P is being calculated. This point of course is fixed and its co-ordi- 
nates constants so far as the above integration is concerned. 

When it is desired to introduce polar co-ordinates explicitly into 
the above and subsequent integrals we know that 

dY' = r'hinO'dr'de'dcl^' 

r^^=.r' + r'2 - 2rr'{cos(9cos^' + sin(?sin^'cos(</) - ^')}. 

Lamb in the above-mentioned article also shows that if 

f + ^ + ffJL = (35) 

dx dy dz 

throughout the fluid, then L, M, N are like P expressible by in- 
tegrals representing potentials. This condition holds if over every 
surface bounding the fluid the vorticity either vanishes or has its 
axis in the tangent plane to the surface. To represent these in- 
tegrals in a form suitable for our present purpose, let 12 represent 
the resultant vorticity at any point and (PQ) the angle between 
any two directed quantities denoted by P and Q. Then the in- 
tegrals are 

L = 



2 



i_ r r r n'cos{n 'x) ^y,\ 

'^^JJj n "" > (36). 

etc. J 

We have now to find the equivalent in polar co-ordinates of 
equations (35) and (36). 

Replacing differentiations with respect to x, y, z by differentia- 
tions with respect to r, 0, </>, and having regard to (32) after differ- 
entiation, we transform (35) into 

i.^(Rr») + _J_;^(S8in^) + -4-,g = ... (37). 
r^ dr rsinOdO rsmd defy ^ 

This result may also be obtained as follows. Suppose L, M, N 

to represent the components in the directions of the Cartesian axes 

of a certain vector quantity. Then under the conditions existing in 

(35) the vector quantity must satisfy the same condition as the 

velocity in an incompressible fluid. But by (32), R, S, T are the 

components of this same vector quantity along the fundamental 

directions at the point r, 0, <f>. Thus they must satisfy the same 

condition as tlie velocity components in those directions in an in 

compressible fluid. Thus (18) should be satisfied when 
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_^ = 0, w = R, r = S, w = T \ 
8t 

and when these substitutions are made we simply get equation 
(37). 

We likewise easily find from (32) and (36), noticing that and 
</> are constants so far as the integrations are concerned, 



T=irf r"'*^"'("'f) ,/v' 

■2^} J J n 



(38). 



If we denote by K the vector quantity whose components are 
L, M, N, or R, S, T, we obtain at once from (38) 

R^ + St; + Tf = Kficos(Kft) 



27r 



"12fi'cos(i212') ^y, 



(39). 



Thus 



f r rK12cos(Kft)c^V = i- 



«/«/•/ 



flfl-cos(12 12^) ^y ^, ^^^y 



The meaning of the sextuple integral is that the value of a 
certain function for the fluid occupying an element rfV at the point 
(r, 0, (j>) is derived from an integration throughout space, and that 
the value of this function is then integrated throughout space. 

When the fluid is incompressible and unlimited, being at rest at 
infinity, and all the vortices are within a finite distance of the origin, 
Lamb shows in his Art. 135 that a certain expression which is 
obviously identical with the right-hand side of (40) is equal to E//), 
where E is the total kinetic energy of the fluid motion. 

It may be seen from the same and previous articles that when 
the. fluid, remaining incompressible, is not unlimited, the right-hand 
side of (40) still equals E/p provided over the boundary or boundaries, 
S, the axis of the resultant vorticity is in the tangent plane to the 
surface, and 



r 


• 


X. fl V 




Kq 


I III n 


-^ 




I' m' 7l' 



(/S = 



(41); 



where (A, /x, v) are the direction cosines of K, (Z, m, n) the direction 
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cosines of the normal, and (/', m\ n) the direction cosines of the 
resultant velocity q. 

Supposing the above conditions satisfied when there exist a series 
of vortex filaments of strengths m^ . . . m^, m^ ... m„ forming closed 
curves, «i, ... Sp, s^ ... «„ in an incompressible fluid of density p, we 
obtain for the energy E of the fluid motion the equation 

= p ^"l7n^\Kcos(Kd8j,)ditj, I ... ... ... (42). 

In the first summation ds'j, is an element of the same vortex 
filament as ds^, and r^ is the distance between them. In the second 
summation every pair of the n closed curves must be combined 
together, so that Jn(n- 1) integrations are included. In the last 
summation K must of course be regarded as a function of the arc ^^ 
of the curve along which the integration is being taken. The value 
of K depends on the combined action of the vorticity in the curve 
over which the integration is being taken and of the vorticity in 
each of the other n-l closed curves. Basset gives in his Art. 92 
formulae which are equivalent to the first of the two expressions given 
above for the energy ; and points out the analogy of the first in- 
tegrals to the co-efficients of self-induction and of tlie second integrals 
to the co-efficients of mutual induction of a series of closed electric 
currents. 

The similarity of the equations of vorticity to those of electro- 
dynamics * may, I think, be brought out very clearly by the follow- 
ing comparison of (39), (40), and (42) with corresponding electro- 
dynamical expressions. We require first to notice that according to 
the table in Basset's p. 88, the following are analogous quantities in 
the two subjects : — 

Vorticity 12 and electric current (g, 
K and electromagnetic momentum til. 

Suppose with Maxwell in his Electricity, Art. 616, that the 
quantity he terms J is zero. Then treating equations (5) f of tliat 

* Cf. Basset's Art. 95. 

t Second Edition. In the first edition change fi into l/fi in these equations and 

in Art.«617. 
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article precisely as we treated our equations (38), we obtain as the 
equation corresponding to (39) 

m(Bcos(iS(E) =/*[!"[ ^^'"^^^^'^ dV' ... (43). 

The meaning of the dashes and of c?V' and r^ is exactly the same 
as in (39). The co-efficient fi is unity in the electromagnetic system 
and denotes the permeability in the electrostatic. In obtaining the 
equations from which (43) is derived Maxwell apparently treats fi 
as constant throughout every portion of space wherein currents are 
actually flowing. 

From (43) we find the following equation answering to (40) 

J I U-'w®cos(SB«)rfv= r rrrj[««:?^««:)^vrfv'... (44), 

the meaning of the sextuple integral being the same as previously. 

Suppose now that the current system consists of currents of 
strengths i^ ... i^, % ... i„ traversing closed linear conductors Sj ... 
8j„ Pg ... 8n oi small cross section, then we may replace 

(B dV by i ds. 

We thus transform (44) into 

i2« ^l J J ^.J^^V) ds, as', j ^. V |- .^ .^ J Jcos(^., dsg) ^^^ ^^^ J 

= P« Fi^ L-i®cos(®(/5^) c/5^ 1 (45). 

The limits of integration and the meanings to be attached to the 
letters are the same as in (42). As /x = 1 on the electromagnetic 
system the second side of (45) is, it will be seen, identical with the 
expression given in equation (14) of Maxwell's article 634 for the 
electrokinetic energy in the field. The first side is also a known 
form for the energy, the integrals of the first series being the co- 
efficients of self-induction, and those of the second series the co- 
efficients of mutual induction of the several circuits. The resem- 
blance of (45) and (42) both in their final form and in their mode 
of derivation is very close. 

I am not aware that any name has yet been assigned to the 
quantity K. It will appear, however, from the previous investiga- 
tions that it is a quantity of considerable importance. It might with 
considerable fitness, as showing its analogy to the corresponding 
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electromagnetic quantity, be termed the vector-potential of vorticity. 
This name would also fit in with the mode of derivation from the 
vorticity, as the following geometrical interpretation of equations 
(38) will show. 

We take the mean vorticity of each element of volume and divide 
it by the distance of the element from the point at which K is to be 
measured. Then starting from some fixed point as origin we draw, 
as in forming the polygon of forces, a series of lines representing 
these quantities in magnitude and in direction parallel to the corre- 
sponding vorticities. The line joining the origin to the end of the 
last line so drawn represents K in magnitude and direction. 

This is practically the same construction as Maxwell gives for 
the corresponding quantity 5 in his Art. 617. I fail, however, to 
see how he obtains the equation 



81 = ;a \ — dxdydz 



in the same article. It seems to me clearly wrong unless all the 
currents are parallel. Employing our previous notation, I believe 
the correct forms of this and of the corresponding hydrodynamical 
equation are 

l& = ^\\\^m)dV' (46), 

K = 2lJ|J^^).;V' (47). 

In each case there exist two other equations determining the 
direction of the vector potential. 

If in the equations (28a) - (30a) for a non-rotating fluid we make 
w, ^ and 7/ vanish, and suppose along the vortex element r, u and v 
constants, we get the case of a circular vortex ring whose centre lies 
in the line ^ = and whose plane is perpendicular to this line. The 
equations show that ^/prsinO must be an absolute constant, and this 
may easily be seen to follow from the equation of continuity. Such 
a vortex ring, as I shall immediately show, can exist outside a sphere 
of which the line ^ = is a diameter. Assuming this I shall first 
consider whether such a ring can exist when the fluid is rotating 
about the line in which the axis of the ring lies. 

As in the similar case in Cartesian co-ordinates, it may easily be 
found that the equations (28) and (29) cannot be satisfied by sup- 
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posing w, ^, 7) to vanish, and u and v to be independent of </>. Sup- 
pose, however, that when w is zero these equations are satisfied by 

^ = 7; = C=Co) • ^ 

Then when w ceases to vanish they will be found satisfied by the 
following, which are consistent values, 

u = u^ v = t\^ w— - (orsin^ \ (x^\ 

supposing p to have the same value as when w vanishes. 

The method of proof is the same in all three of equations (28) - 
(30), so it will be sufficient to exemplify it by treating (28) only. 

Substitute the values (49) in (28), and then noticing that when 
0) is zero (49) satisfies the equation, we are left under the necessity 
of proving 



8 



+ -!^1 -77^. sin^(woSin^ + VqCOs^) \, 
sinc7 \(W J 



Substituting for ^ from (18), and then for — - from (17), we 

p ot ot 

have to prove 

[cos(9 1- - ?1!L I {r-u^) + _^ I — . sin(9(wosin^ + t7ocos^) |. 
dr r ad A sm^ \dd J 



= 0) 



Carrying out the differentiations this will be found to be an 
identity. Similarly it may be shown that (49) satisfies (29) and (30). 
In the last case indeed the fact is obvious on inspection. This in- 
dicates that so far as the ring itself is concerned the equations may 
all be satisfied if the fluid in the ring possess, in addition to the 
velocity it would possess in a non-rotating fluid, a rotation with 
uniform angular velocity w about the axis in an opposite direction 
to that in which the fluid and axes are rotating. Each material 
cross section of the ring in fact would move in a plane fixed in 
space. 

The difficulties that occur in the motion of the fluid surrounding 
the ring are precisely the same as in the corresponding case in Car- 
tesian co-ordinates. We may, as there, by bringing in the considera- 
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tion of the existence of viscosity in all actual fluids, arrive at a 
similar conclusion as to the probable character of the phenomena in 
nature. 

In order to prove that the motion we have just investigated can 
apply to the case of a vortex ring when the fluid is bounded inter- 
nally by a spherical surface, a diameter of which forms the axis of 
rotation and coincides with the axis of the vortex, we still require 
to show that all the requisite conditions for a vortex ring in this 
position can be satisfied when the fluid and axes are not rotating. 
All the data necessary for this already exist. 

It has been shown by Mr T. C. Lewis* that if a vortex ring of 
strength m exist in an incompressible fluid outside a spherical 
boundary of radius a, with its core at a distance ^ from the centre, 
the boundary conditions are satisfied by the velocities deducible from 

the true ring and from an image ring of strength - m Jfjf' inside the 
sphere with its core at a distance /' = a7/ froni the centre. Corre- 
sponding points of the tAvo rings are, it will be noticed, in the same 
relative position as an electrified point and its image with respect to 
a spherical conductor. 

Again it has been shown by Mr W. M. Hicks f that a " source" 
of strength m at a distance y*> a from the centre O of a spherical 
surface of radius a has an image inside the sphere, consisting of a 
"source" of strength nia/J'a.t the point Q which is the inverse of the 
true source P, and of a line " sink " of strength ni/a per unit length 
extending from Q to O. 

Now by a source of strengtli m at P we simply mean that there 
exists a term - m/r^ in the velocity potential, where r^ is the distance 
from P. But from (34) if a vortex ring, whose element of cross 
section is da- and element of core ds, exist in a compressible fluid it 
contributes to P, which there denotes tlie velocity potential, the 
expression 

47r J \J J p' 86 r, I 

The double integral inside the bracket extends over o-, and the 
integration with respect to ds is taken round the core. 

Assuming, as usual, that the radius of tlie cross section is very 
small compared to by the radius of the ring core, we may regard r^ as 

* Quarterly Journal, Vol. xvi., p. 338, or Basset's Art. 311, Vol. ii. 
t Phil. Trans. 1880, Part ii., p, 455, or Basset's Art. 52, Vol. i. 
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constant over the cross section and equal to the distance of the 
external point from the core. We may then replace the above 
expression by 

47rJ Vi 

JJ P 8t 
I have shown on a previous occasion* that when p is uniform 
over the cross section, b being very large compared to the radius of 
the section, we may take 

8t 

Whether p vary or not over the cross section, provided & be the 
same for every point on the core, a ring of this kind is equivalent to 
a ring source of uniform strength 6-/47r per unit length of core. 

Applying Mr Hicks' result, denoting by a the radius of the 
sphere and by / the distance of the ring core from the centre, we find 
for the image a ring source at a distance or// from the centre, whose 
strength per unit length of core is 

«^/^ ^ = Z_ ^ 
y or 47r a 47r' 

and a sink spread over the surface of a cone which extends from the 
image source to the centre. The strength of the annulus of the sink 
by two planes perpendicular to the axis of the cone intercepting unit 
length on the generators is everywhere 

27rb& h . 



cr. 



47ra 2a 

The strength of the cone sink per unit of surface would thus 
become infinite at its vertex, but that only means that there is a 
point sink of finite strength at the centre of the sphere. Combining 
the systems of images for the vorticity and for the compressibility, 
we clearly possess a solution of the problem whose solubility we had 
to establish. 

While the preceding investigation as to the image system for the 
compressibility possesses the advantage of bringing out more clearly 
the physical side of the question than is possible in the analytical 
method, it would not in general lead very readily to the actual de- 

* Proceedings^ Vol. vi., p. 65. 
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termination of the fluid velocities. When the angle 2 a of the cone 
subtended by the ring at the centre of the sphere is small, we could 
obtain approximate values for the velocity at points at a considerable 
distance from the ring by regarding it and the image ring as point 
sources and the image sink as a line sink. 

In general recourse had better be had to the analytical expression 
for the velocity potential given below. Using a, 5, a, and f in the 
same sense as above, we have h =ysintt. 

Denoting by Q„(^) the nth zonal harmonic whose pole, ^ = 0, 
lies on the diameter of the sphere which is the axis of the ring, I find 
for the complete value of the velocity potential at the point (r,^) 
referred to the centre of the sphere as origin 

where 



^{j-^Q.(a)Q.W+... 



n 



-^Q»(«)Q«(^)+-} - (50), 



p,= - 



n + 1 /»- V+i 
-^{l + f Qi(a)Qa(e)+ ... +^Q„(a)Q„(e)+ ... } 



. ;^(5i). 

= - -^f A + <Q,(a)Q,(e) + ... + ^Q„(a)Q„(^) + ..) 

The first paft of the velocity potential, P^ may be regarded as 
answering to the image system. The second part P2 is the velocity 
potential due to the ring itself ; and the first or the second value for 
it is to be taken according as the point considered is nearer the 
centre of the sphere or is more remote than the ring core. 

In calculating Pa the method suggested in Thomson & Tait's 
Natural Philosophy ^ Art. 546, may be readily applied. The deter- 
mination of Pi is based on the vanishing of the normal component 
of the fluid velocity over the spherical surface. It presents no 
difficulty. 

Over the spherical surface the velocity is everywhere of course 
along the tangent which lies in the plane through the axis of the 
vortex ring. At an angular distance from the axis its magnitude 
is 

[1 /c?P, . c/Po\l o-sinasin6' r3«>^ / \r\ >/ \ 

+ ^ (y )"Q-.(«)Q»'M + • • • } •• • (52), 
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where Q,/(/a) = — Q/,(a''), denoting cos^ by /x, and replacing the nota- 

tion Q„(^) by Q„(/x). The surface velocity thus vanishes, as is 
obvious from symmetry, where the axis of the ring cuts the spherical 
surface. 

The expressions for the velocity potential in the form given above 
do not converge rapidly unless the distance of the ring from the 
centre of the sphere be considerable compared to the radius. 



On a problem in permutations. 
By R. E. Allardice, M.A. 

The problem to be considered may be stated as follows ; — How 
many necklaces may be formed with p pearls, r rubies, and d 
diamonds ?* 

The peculiarity of this problem is that a general solution cannot 
be given in terms of p, r, and d alone. The form of the solution 
depends on the nature of the numbers p, r, and d ; and it is neces- 
sary in solving the problem to consider whether or not these num- 
bers have a common measure, and how many of them are odd and 
how many even. All possible cases of the problem are not discussed 
in this paper ; but enough of them are considerod to illustrate the 
variety of forms that the solution may assume. 

If we put p-{-r + d = n, the number of possible arrangements of 
the n stones in a line is n\/p\r\d\ Hence the question is, how many 
of these arrangements will give the same necklace ; or, conversely, 
if we take any one form of the necklace, how many different 
arrangements of the stones we can get from it by breaking it at 
different parts and stretching it out straight. It is obvious that if 
the n stones had been all different, the answer to the second of these 
questions would have been 2n ; in other words, with n stones all 
different, we may form nl/2n necklaces. The further question then 
naturally arises, In what cases, if any, are the 2n arrangements of 
the stones obtained from each form of the necklace all different when 
the stones are not all different ? Now these 2n arrangements com- 
prise the n that are obtained by a cyclical interchange of the stones, 
one at a time, together with the n that are obtained by exactly re- 
versing each of these n arrangements. 

* This problem was suggested to me by Professor Chrystal. 
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The 2n arrangements may fail to be all different for one of the 
three following reasons : — 

(1). Because an arrangement is not altered when it is reversed. 
In this case the arrangement must be symmetrical. 

(2). Because an an*angement is reproduced after a cyclical inter- 
change of a certain number of the stones. It may easily be seen 
that in this case the arrangement considered must be resolvable into 
a number of identical groups. 

(3). Because an arrangement is reproduced when a number of 
the stones have been interchanged cyclically and the resulting ar- 
rangement reversed. In order that this may happen, the arrange- 
ment considered must consist of two symmetrical groups. These 
three cases may be conveniently referred to as the case of a single 
symmetrical group, the case of identical groups, and the case of two 
symmetrical groups, respectively. 

A single symmetrical group will occur if not more than one of 
the numbers p, r, d, is odd ; and not otherwise. 

Identical groups will occur if p, r, and d have a common measure ; 
and not otherwise. 

The case of the two symmetrical groups will occur if not more 
than two or the numbers jo, r, (/, are odd ; and not otherwise. 

First Case. 

Hence if p, r, and d are all odd and have no common measure, 
the number of necklaces that may be formed is n\/2n.p\rld\f but this 
formula will not hold good in any other case. [It may be pointed 
out that any two of the numbers may have a common measure and 
the formula will still hold ; and that a similar formula will apply to 
the case of stones of any number of different kinds, provided there 
be an odd number of each of three of the kinds.] 

Second Case. 

Suppose now that two of the numbers, p, r, d, are odd ; but that 
these numbers have still no factor in common. 

Let;? = 27r+1, r=2p+l, d=28. 

We may now have two symmetrical groups, a pearl occurring at 
the centre of one and a ruby at the centre of the other. It should be 
noted further that a cyclical interchange in the case of two symmet- 
rical groups produces two symmetrical groups ; and that by means of 
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a number of cyclical interchanges one of the groups may be reduced to 
a single stone. Hence, if we put v = tt + p + 8, the number of neck- 
laces involving two symmetrical groups will be vI/ttI/jIS!. From each 
of these may be obtained n linear arrangements of the stones ; and 
thus there are altogether n\/p\r\d\ - n, vl/irlplSl arrangements that do 
not involve symmetrical groups, giving (n\/p\r\d\ - n,v\/'jr\p\8\)/2n 
necklaces. Hence the whole number of necklaces is 

1 I w! _ n.vl I *'• _ (^"^) . *'• 

2^ I plrldl 7r!/^!j tH^" 2.p\r\d\' 2^d^! 

It may be well to illustrate this formula by actually writing out all 
the possible arrangements in a particular case. 

Put p=l, r = 3, d=2; then 7r = 0, p=l, 8 = 1; ?i = 6, v=2. 

The formula gives '- — + = 6. 

^ 2.1!3!2! 2.0!1!1! 

The arrangements are the following, the two in the first column 
containing, and the others not containing, symmetrical groups : — 

prdddr prrddd pddrrd 
pdrdrd pddrdr prddrd 

Third Case. 

Next let only one of the numbers p, r, rf, be odd ; and suppose 
that these three numbers have still no common measure. 

Let jo = 27r+l, r = 2p, c?=28, n=^ + r + (/, v = 7r + /j + 8. We 
have now the case of a single symmetrical group, and also the case 
of two symmetrical groups to consider. 

In the latter case one of the groups will contain an odd number of 
stones and will have a pearl as the centre one, while the other group 
will contain an even number. It may easily be seen that any such 
arrangement may be reduced by a number of cyclical interchanges 
to a single symmetrical group. 

The number of necklaces involving symmetrical groups is i'I/ttI/jIS!, 
giving n.v!/7r!/3!8! different linear arrangements. Hence it may easily 
be seen that the whole number of different necklaces is 

(n-l)!/2.jt>!rW! + i/!/2.7r!/o!a! 

This formula is exactly the same as that obtained in the last 
case ; but ir and v have here slightly different meanings. 

As an example, put jo = l, r = 2, c?=2; so that 7r = 0, p=l, 8 = 1 ; 
71 = 5, v = 2. 
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The formula give8 4!/2.1!2!2! + 2!/2.0!l!l! = 4. 
The arrangements are the following, the two in the first column 
containing, and the other two not containing, symmetrical groups : — 

rdpdr rrpdd 

drprd drdrp 

Fourth Case. 

Suppose in the next case, that p, r, and d are all even. 

This case is an exception to the assumption made hitherto that 
p, r, and d have no common measure. I shall consider only the case 
in which the other factors, when 2 is divided out, are all odd num- 
bers, and have no common measure, that is, 

p = 2'jr, r = 2p, d=28; where tt, p and 8 are odd numbers and 
have no common measure. 

We may now have (1) a single symmetrical group ; (2) two 
symmetrical groups ; (3) identical groups. 

Under (2) we have to consider two cases, namely, (a) that in 
which each of the symmetrical groups contains an odd number of 
stones, the central stones of the two groups being necessarily of the 
same kind ; (fi) that in which each of the symmetrical groups con- 
tain an even number of stones, an arrangement which is reducible to 
a single symmetrical group by means of cyclical interchanges. 

(1). The number of necklaces involving single symmetrical groups 
is v!/2.7r!p!8!, giving w.v!/2.7r!/o!S! permutations. 

(2). The number of necklaces involving two symmetrical groups 
of type (a) is 2{(i/- l)!/2.(7r - l)!p!8!}, the three terms corresponding 
to the three cases where the middle stones are 2 pearls, 2 rubies and 
2 diamonds. The number of permutations that these give is 

w.2{(v-- l)!/2.(7r- l)!/o!8!} = n.(v- l)!27r/2.7r!/)!a! = 7i.v!/2.7r!p!6! 

(3). The number of permutations involving identical groups is 
v!/7r!/o!8! giving (v- l)!/2.7r!/o!6! necklaces. 

Hence the whole number of permutations not involving any of 

these three cases is 

n\ n.v\ v! 

and the whole number of necklaces is 

1 / n\ n.v\ v\ { v\ (i/-l)! 

2^ I ^Mdl Jr\p\8\ 7r\p\8\ J "^ wlplSl "^ 2.7r\p\8\ 

^(n-J^ v!_ (.-1)! 

2.p\rld\ 2.7r\p\8\ i.7r\pl8\ 
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As an example, put p«2, r=a2, d=^2; so that 7r=»l, p=l, 
S = l, w = 6, v = 3. 

The formula gives 51/2.2.2.2 + 3!/2 + 2!/4 = 1 1. 

The arrangements are the following, those in (1) containing a 
single symmetrical group, those in (2) two symmetrical groups, 
those in (3) identical groups, and those in (4) having none of these 
peculiarities : — 



(1) dprrpd 
pdrrdp 
prddrp, 



(2) pdrprd\ 
drpdpr \ 
rpdrdpj 

(3) pdrpdr 



(4) prddpr\ 
rdpprd I 
pdrrpdi 
pprrddJ 



Fifth Cask. 



Suppose next that 77, r, and d have a factor in common. In a 
full discussion of the problem a number of particular cases would 
require to be considered ; but I shall limit myself to that in which 
/>, r, and d are all odd numbers and their G.C.M. is a prime number. 
In other words, I assume p = Att, r = A./3, d= A.8, where A, tt, p, and 
8 are all odd numbers, A. is a prime, and tt, /a, and 8 have no factor 
in common. 

A permutation may consist of A, identical groups, each containing 
V = TT + p + 5 terms. 

There will be v!/7r!/o!8! such permutations, giving (v!/7r!/3!8!)/2v 
necklaces. 

Hence the whole number of necklaces is 

_l_(_n! 1^1 (v->l)! ^ (n->l)! (n-v).(v-l)! 



,p\r\d\ 
2.p\r\d\ 



A-1 



2.7r!p!5! 



As an example, put jt> = r = c?=3; so that 7r = /3 = 8=l, n = 9, 
V = 3. The formula gives 8!/2.3!3!3! + 2/3 = 94. 

Only one of the necklaces contains identical groups, namely, 
prdprdprd. 

Without writing out all the 94 arrangements, we may verify that 
this is the correct number by enumerating them in the following 
manner : — 

We may set down the three pearls in a ring, and count the num- 
ber of ways in which the other six stones may be arranged relatively 
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to them. The number of stones to be put in each of the three spaces 
between the pairs of pearls is given in the following table, with the 
total number of corresponding necklaces in each case : — 

6 gives 10 necklaces, 
20 
20 
10 
10 
20 
4 

94 
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2 
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Sixth Meeting, 2\th April 1890. 



R. E. Allardice, Esq., M.A., Vice-President, in the Chair, 



On a hydromeclianical theorem. 
By Dr A. 0. Elliott. 

Giffard's injector appeared more than thirty years ago. The 
first serious attempt to explain its action on dynamical principles 
was made by the late William Froude at the Oxford Meeting of the 
British Association in 1860. The history of mechanical science is 
almost everywhere deeply marked by Rankine ; and it seems, just 
as it ought to be, that he should be found to have contributed not a 
little to the literature of this particular subject in a paper presented 
to the Royal Society of London in 1870. As serving to show how 
far the problem is still interesting, even from a high standpoint, 
attention may be directed to the exceedingly curious procedure of 
Professor Greenhill, where he deals cursorily with the matter at the 
page numbered 448 of his article on Hydromechanics in the Ency- 
clopcedia Britannica. 

When first announced, the statement that the particles of a mere 
steam jet could, by the agency of this somewhat simple apparatus, 
force for themselves, in addition to a considerable quantity of more 
or less cold feed water, re-entrance into the identical boiler from 
whence they had escaped, seemed to involve an impossibility. But 
the mystery of that aforetime paradox would have been as nothing 
had it then been farther known what is now f amiliar— namely, that 
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the steam for the steam jet might be in some cases taken from the 
exhaust-pipe, instead of the boiler, with equally effective results. 

In dealing with the theory of Giffard's injector, the author found 
that it was convenient to employ the principle of impulse and 
momentum ; by which is meant simply a statement of Newton's 
second law in the form that the time-integral of the forces acting on 
the masses concerned, and reckoned in a specified direction, is equal 
to the resulting gain or destruction of momentum reckoned in the 
same direction. 

In the ordinary symbols but merely altered in respect that y and 
V are regarded as component vectors, Newton's second law is repre- 
sented by /. dv 

and hence 

= mdv. 



fa=j 



The idea is applied in hydromechanics by considering a fluid- 
filled closed surface ; and then the sum of the time-integrals of 
the external forces acting on or through the bounding surface 
is equal to the generation of momentum, all reckoned in one and 
the same direction. It is an instance of Mr Froude's remarkable 
insight that, in the paper referred to above, his explanation of 
Giffard's injector was based, so far as it went, on these principles. 
And it is worthy of remark, too, that Professor Greenhill, in the 
article already mentioned, by the help of Green's theorem, deduces 
Euler's equations directly from the fundamental principle. 

But first we must obtain the equation of energy. Suppose a 
fluid to move steadily, irrotationally, and frictionlessly along a tube 
of stream lines. Let Ai, jt?„ Vi and pi be respectively the area, the 
pressure, the velocity, and the density at one section ; and let similar 
symbols with the subscript 2 represent the same quantities for 
another and second section. Farther, let q be the volume corre- 
sponding to the mass-flow Q and density p. Then neglecting for our 
purpose differences of level, and equating the work done by the 
pressures on the end sections per unit time to the increase of energy, 
we have 

where the pressures are measured in terms of the weight of unit- 
mass. The equation of continuity gives 
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AiVipi = A^V2pi = Q; (2) 

and qp = Q, ... ... ... (3) 

\pdq = pq- \qdp, ... ... (4) 

\ pdq = A^V2P2-A^v^Pi-Q\ -^. ... (5) 

Ji J^ P 

,\ (1) becomes 

= f-f-+r±; (6) 

^g '^g }^ p 



or 



^+r^ = const = H (7) 

^g } P ^ ^ 

For the case of a liquid p is practically constant ; and the equa- 
tion of energy becomes 

— + .^ = 11, ... ... ... (o) 

2^ P 

where H may be called the total head. 
From (8) 

;? = ^|H-_| (9) 

The equation of continuity gives 

Avp = Ai^i/o = A2V2P = Q. 

^A^-hM- ■ ■ ^''^ 

Multiplying by dA and integrating the impulse of the sides of 
the tube per second 

.JVa=4h(a,-a.)*^(1-1-)] 

= 2/oH(A2- Ai)+jt?iAi-/?2A2. ... ... ... (11) 

Now, considering the impulse on the areas Ai, A2 it appears that, 
to obtain the total impulse, there falls to be added to the last 
expression 

;?iAi-;?2A2. 
Hence the total impulse per second, say 

J = 2{piAi-j03A2-pH(Ai-A,)} (12) 
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= 2(p,-p,)-^ (14) 

(13) may be obtained ahnost directly ; for since impulse is equal 
to momentum generated 

9 

as before. 

Take for an illustration of (14) the case of a liquid escaping from 
a large reservoir by a bell-mouthed tube into a space where the 
pressure is zero. Then Ai = go , jo^ = 0» Therefore 

and not simply PiA^ as might be hastily and erroneously concluded. 
Secondly, suppose the fluid to be a gas and that the relation of 
pressure to density can be represented by 

p = kp^ ... ... ... (16) 

where k and y are constants. 
Reverting to (7) we have 

'dp y p 



P 



P 7-i P 

^^ /v ■ * ••• ••• ■•• \ J 

P 



where A. = — ^- — 
7-1 
The gain of momentum is as before 

9 
and the equation of continuity gives 

Q = Avp = const. 

, *. tl = A.2P2 — "^iPl * ••• '•* "•• \/ 

9 9 

Substituting from (17) in (6), and by the help of that result and 
(2) eliminating v^ and Vg from (18), 

J=.2\II^-^]^'PliA^ (19) 

(14) and (19) are of course analogous ; and the first may be 
obtained from the second by putting /5 = const. The author has 



73 

found these expressions to possess some interest ; and it was princi- 
pally with the object of drawing attention to them that the present 
paper was' entitled as above. 

Let steam escape from a reservoir at pressure pi into a space at 
pressure p^ ; and suppose the law 

p = kp^ 
to hold. Then from (6) 

!!L=r^ = x(j^-^\; (20) 

or, v= J2f\ l?l-h. 

V Pi Pa 

Hence the mass-flow 



Q = Aap, V^i^'A. jh.-I!l (21) 

V Pi p'i 

Q is maximum when for a fixed value of pi^ p^ has such a value 
that 



That is when 



Therefore 



^p\-p2P2 a maximum. 
Pi 



J»2 _ 2 ;?, 



Pa y + 1 Pi 



(22) 



/y-1 lih 



It is found by experiment that this maximum value of Q occurs ; 
but that on continuously diminishing p^ below the point correspond- 
ing to the maximum hardly any appreciable change takes place. 

3 

In fact where j92< — P\y ^.n expression of the form (23) holds ; and 

o 

consequently 

Q = mA2^/>i/3, (24) 

where m is a constant not much differing from 
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Accepting the guidance of these facts, and reverting to (19), we 
find when Aj = oo that 

J = 2klPl-£l\A^, (25) 

This expression will be maximum for a given value of ^ when 

Pi 

^P.-P. 
Pi 

is maximum. That is when 

Pt 7 Pi 
Hence the maximum value of J is 



2 
Pi 

= 2^A2P2 
Pi 

^2(~)y=^p,A, (26) 

= 2??;PiA2 ... ... ... (27) 

where n is put for / ^ \7iiT. We shall take this as a probable value 

of J for the steady state. 

Fig. 55 represents diagramatically the essential parts of an in- 
jector. S is the steam pipe ; W, the water inlet ; O, the overflow ; 
and B the pipe leading through the check valve to the boiler. Let 
Pi be the boiler pressure ; A, the area of the steam nozzle ; and A3 the 
area of the water-steam throat. The author intends on some future 
occasion to make an experimental determination of the pressure in 
the chamber where the steam and water mix. In the meantime, if 
the feed water be not too hot, we may assume that the pressure is 
but little greater than zero. 

Let us suppose also that the water and steam pipes are large, 
and that the passages and nozzles are well and gently tapered. 
Then in applying (14), Ai may be taken as =00. Farther, let v^ 
represent the velocity of the water-steam jet, p^ its pressure, and p 
the density ; and put p^ for the atmospheric pressure (diminished by 
a quantity corresponding to the lift of the injector, if any). Then 
taking the sum of the impulses and equating to the gain of mo- 
mentum, there results 
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v' 



2/AP1A2 + 2p^Acosa -p^As -/i/j^Ag = A3V3P , Vs\g, ... (28) 

where /^ is a friction and eddy resistance factor, A is the sectional 
area of the water passage projected on a plane at right angles to the 
stream lines, and a is the angle which a tangent to the stream lines 
at the same place makes with the axis of the injector. Now, from 
the equation of energy, 

Ps Pi vl ,rvl, 

p p 2g 2g 

^^ ^=^-/,ji (29) 

P P '^9 

where /o is a friction factor, and/j is written for 1 -/q. 

Hence putting Acosa = sAj 

^ _ 1^9 l 2(np^ + 8p,)A^-p^A^ .oAx 

"''Ja^ p{i+A)As - ^^^> 

The mass of steam used per second is say 

S = mAj JpiPi ; 

and the mass of water passed on to the boiler per second is 

A3V3/). 
Therefore the mass of water injected per second, say 

W = A^VsP - mA^ Jpip^ (31) 

Hence the ratio of the mass of water injected to the mass of 
steam used is 

Y^_ A,v,p_ ^^ ^3^^ 

S mAj Jp^pi 

The efficiency of the apparatus will be maximum when 

mAo JpiPi 
is maximum ; that is when 

A, 

— -Vo 

A/ 
is maximum. Writing r for the ratio A3/A2, it appears from (30) 
that there is to be made maximum 

2(npi + »jt?„)r-/?r*. 

Hence for best working 

^^ np, + sp, (33) 

Pi 
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This value of r makes 



Pi= . 2f ^^ (^^) 

If y*! were small enough, then this condition could not be fulfilled 
without 

P^<Pa 

which is inadmissible with an open overflow. 

On the other hand, if /i were large enough, this value for v could 
only be consistent with pressure above atmospheric in the overflow 
chamber, which, under the same conditions as formerly, is equally 
inadmissible. 

It appears, therefore, that the best adjustment will tend towards 
the production of pressure in the overflow chamber above or below 
the atmosphere according as 

-^ /?! IS > or <;?„. 

With an open overflow we must have 

y 2 K~ ^^ • ••• ••• ••• V*^^/ 

% P 

Therefore from (28) putting 773 =^9^, 
or ^ ^ 2(71/?, + sp^) 



'^ (Pl -Pa) +P. 



(37) 



(38) 



It appears also from (28) that in any case 

r must be < 1(!'^1±^ 

Pi 

If Tj and T, be respectively the temperatures of the steam and 

of the feed water in degrees Fahr., then in order that there might 

possibly just be condensation at the throat. 



^_mV;,.p j912 + -3(T.-32)L „^„. 

but the mutual action is not actually so perfect, and it is well to 
allow 
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8 



^m>iP^{974 + 3(T^->32)} 



(josa. ... (39) 



Taking 7= 10/9 

w = l — J7^=:'4 nearly. 

Suppose, further, for example, that the boiler pressure is 160 lbs. 
per sq. in. absolute ; that the diameter of the water-steam throat is 
•25 inch ; and assume the following values : — 

T,= 60'*F.; a = 13^/=-5; 

/2= -9 ; m= 36, allowing for contraction. 

Then r = 42 ; and diameter of steam nozzle 

= '^^ ='38 inch. 



V-42 

f m ^ 

W = A,( j2gp(p,-p^)//2-- J,,,p,j ... (40) 

= 3*0 lbs. per second. 

For the exhaust steam injector we have merely to substitute p^ 
forpi where it is associated with n, and in (31) where it appears 
under the radical. Thus 

'^Pa(n + s) 

^ = 2+/i/ v^ (41) 



8 



J2gp{l50' -T,} 



(42) 



and . , .- — : r_ m 



W = Aal J2gp{p, - pM - - ^PaPah . . . (43) 



On Bazikine*8 Formula for Earth Pressure. 
By Dr A. C. Elliott. 

In a short course of lectures on " Railway Practice," the author 
was recently called upon to deal with the mechanical principles 
involved in the design of retaining walls. What has come to be 
known as Bankine's method had to be explained, at all events, in its 
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practical application. But the time at the author's disposal did not 
admit of the general consideration of the theory of stress by which 
Rankine in characteristic fashion leads up (in his Applied Mechanics) 
to the particular problem under discussion. The author had there- 
fore to choose between omitting a demonstration or making a short 
cut, which at the same time should be of a rigorous nature. 

The following process, which in the circumstances was the 
one employed, may be regarded, from a mathematical point of view, 
as a direct solution of the problem : — Given the angle between 
two conjugate stresses to determine their ratio when the maximum 
obliquity on any section has a given value. In the actual problem one 
conjugate stress is produced at any (a certain) point by the weight of 
the superincumbent column of earth ; the other conjugate stress is pro- 
duced by the reaction of the wall ; and the maximum obliquity when 
the wall is exerting its least possible resistance is the angle of repose 
for the earthy material in question. The plane of rupture can of 
course be found, and it will be defined as that plane on which the 
obliquity of the stress has the given maximum value. (Fig. 56). 

Let ABC represent a section of a triangular prism of earth whose 
length-axis is horizontal and parallel to the face of the wall ; and 
suppose moreover that the prism is of unit length. Also let AB^a 
be vertical and parallel to the face of the wall ; and let BO be 
parallel to the upper surface of the earth. Let p, be the intensity 
of the vertical thrust due to the weight of the earth ; and let pj^ re- 
present the conjugate thrust on AB. When, following Rankine, the 
friction of the wall is neglected the direction of pj^ will be parallel to 
BO. Farther, let R be the total intensity of the stress on AC and 
a the obliquity ; P the angle which BC makes with the horizontal, 
and the angle BAO. 

To find BC in terms of a we have 

BC sin^ 

~a~'^cos(^ + ^)* 

Hence resolving vertically and using this value of BC 

Rsin(^ + a)=;7,a — ^^--^^p^asinp; ... (1) 

cos(^ + p) ^ 

and resolving horizontally 

Rcos{0 + a)=pf^acosP (2) 
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putting p^Ip, = S. 

Now for a given value of S to find the plane AC for which a is 

maximum, put 3-5 = 0. Writing a explicitly 

du 

^ + a = tan~^fI._^^^?^|--^-tanM. ... (4) 
•'• 1= ^ ^ 



1 + |1 ^M -tan^r^cos^ 



^ cos(g + j8)cosg + sinfein(^ + ^ ) . 

co&\e+ 13) 



or ^^ 1 1 



1 + (1 ^^ .tan/?y'S^^«^(^ + ^)' - (') 

I S cos^cos((9 + P) ^ J 

Equation (5) determines 6 corresponding to a a maximum. 
When rupture is about to take place the maximum value of a is 
equal to the angle of repose, say <^. To find the value of S corre- 
sponding to this limiting state, we have therefore to eliminate 6 
between the two equations 

and 

Scos»(^ + 1^) n ^ff T^ • (7) 

Is cosy8cos(^ + i«) ^/ 

Using (6), (7) may be written 

ScosX^ + ^) = _J^^_^^ = cos'(^ + <^); ... (8) 

or 7Scos(^ + )8) = cos(^ + <^) (9) 

Hence the two equations between which d is to be eliminated 
may be written 

^Scos)8cos(^ + <^) 
and JScos{0 + P) = cos(0 + €t>) (U) 
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sin^ = sin(^ + </>-</>)> 
expanding and transforming, there results 



tan(^ + <^)fJ^?^^-l)=J_!15t + tani8. ... (12) 



Next putting (11) in the form 



7Scos(^ + <^ - <^ - ^) = cos(^ + ^), 
expanding and arranging, there results 

tan(^ + <^) ^ . . ,^ m {-4 " ^^<^ - ^)} • - <^^) 

sin(</» -p)y J^ J 

Combining (12) and (13) 

sin(<^-MvS ^ ^Vl^scos/^ / ^scos^ 

+ tanj8. ... ... (14) 

Multiplying up and reducing we have finally 

S-2jS-22?^+l = (15) 

COS(/> 



JB'=^± /«^_i (16) 

cos<^ V cos*^</> 

cosl3± Jcos^P - cos'(l> 

And interpreting the signs according to the principle of least 
resistance, we have when the yielding is about to take place by 
failure of the wall 

q^ cosjg- Jcos^P - cos^<j> .,gv 

The other set of signs represents the case of yielding by crushing 
in of the wall, 

[See also Elliott, Proc, R.S.E., Jan. 1887.] 
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Seventh Meeting, 9th May, 1890. 



R. E. Allardice, Esq., M.A., Vice-President, in the Chair. 



On a certain expression for a spherical harmonic, with 

some extensions. 

By John Dougall, M.A. 

The object of this paper is to show how the leading properties of 
Spherical Harmonic Functions may be readily deduced by employ- 
ing as a typical harmonic a certain simple algebraical expression, 
which obviously satisfies Laplace's equation ; and to extend a similar 
method to the case of any number of variables. 

Some well-known expressions of Spherical Harmonics by means 
of Definite Integrals are readily arrived at by this method. 

The most important properties of Spherical Harmonics, perhaps, 
are those connected with the integration of the product of two har- 
monics over the surface of a sphere, and it will be seen that this 
integral takes a somewhat remarkable form when the harmonics are 
of the type I have referred to. 

Consider the function of x, y, z 

u = {ax + by -f- czY 

We have — j -\ -, + - = {a- + Ir + c^){ax + by -\- cc)"~-. 

d'.itr dy'' dz- 

Hence if a^ + Ir + c' = 0, %c is a harmonic of degree n. 



Let a =/+ if\ b = g -V iy\ c = li ■\- ili \ 1= ^ - 1 
The condition a- + h^ + c'- = gives 

/» g^ ^* ; /', g\ ^*' being real. 

Hence by transformation of rectangular axes, we may reduce u to 
the form (x + ty)", or, in polar co-ordinates, r"sin"6^(cosn</) + tsinn</)) ; 
so that the real and imaginary parts of u are two conjugate sectorial 
harmonics. In what follows, then, we consider the general harmonic 
as made up of a sum of sectorial harmonics. 
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Observe that, by a well known proposition, 

(1) i^'*'^\ax + by + cz)~^'^ is another harmonic of degree w, where 

Now let us consider the integral of the product of two harmonics 
of the type u, taken over the surface of a sphere of radius p. 

If the degrees of the two harmonics are positive, but different, 
Green's Theorem shows at once that that integral is zero. 

Let then u„ = (a^x + h^y + Cizy 

v„ = {a^ + b^y + c^Y 

n being a positive integer. 
By Green's Theorem 



\\ 



— u^v^dS 



J J j\dx dx dy dy dz dz / 



and n being a positive integer, we may take the latter integral 
through the whole volume of the sphere. 

But the function to be integrated through the volume being 
homogeneous of degree 2w - 2, we have 
its volume integral = the product of its surface integral into 



I 



— . dr, 

p'"-' P' 



This last integral = - 



2n+l 

Hence w„v„c?*S'=— ^—-(aia^ + 6162 + 0^02) x \u^-iVn-idS 
J J 2n+l J J 

Continuing this method of reduction, we get finally 

(^) 1 h"^"'^^= (2n?i)(2n:i/...3 -^-^^">'"' "^ '^'^ + ''"'^>"- 
And we may also write for reference 

(3) \unVf^dS=0 when, w, n are different. 

In order to apply this result at once, let us find an expression for 
a Zonal Harmonic. We see that (tcosa.a5+ tsina.y + s;)" is a case of 

form, and .'. I (tcos.aa; + tsinay + zyda is a harmonic. 



our 



By putting x= Jc' + y'^cos(t> 
t/= ^x- + 2/^sin</> 



i 
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r2ir 

this becomes I (t Jot? + y'cosa - <^ + zyda 

Jo 

I2ir , 

This function is symmetrical about O2;, and considering aj, y, 2 as 
co-ordinates of a point on a sphere of radius unity, we see that when 
x = y = z~lf the value of the integral = 27r. 
Hence the zonal harmonic with axis Oz, 

1 r27r 

say P«(^) = ^ {tJx^ + y^cosP + zydp. 

J 

Now our result (2) shows that in order to find the surface integral 
of the product of this and another harmonic (ax + by + cz^y we may 
substitute a, b, c for x, y, z, respectively in the expression for ?„(«) 
and multiply the result by a numerical coefficient which we may 
call k. 



We thus get — j '"(, ^^2 ^ j2^jQg^ ^ ^Y^p 





k . r2T 



2^ Jo 



^^ .c« 



2n+l 



Now c** is the value of {ax + by + cz)** at the pole of ?„(«) ; and by 
adding together forms like {ax + by + cz)*" we may get any harmonic. 
Hence the well known important result which we may write 



(4) Jl^^^S=^^^.Y-„. 



Now, considering further the form (ta;cosa+ t2/sina + »)", we see 
by Fourier's Theorem that if we expand it in cosines and sines of 



1 



r2ir 



multiples of a, the coefficient of cossa is — (txcosa + tysina + 2;)" 

TT A 







cos«ac?a. This is obviously a harmonic, and by putting 



2/= Jx2 + yhin<l> 

and reducing, we see that it is a tesseral harmonic of type s. 

The other of the same type is got by writing sinsa instead of 
cos^a in the integral, 
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Now consider the surface integral (r= I) of the product of the 

f r2ir 

above harmonic under the with another (lxcosP -{- Li/sin^ + zy 

J Jo 

coss'/Sdp. 

If 71, 71 are different, (3) shows that the result is zero. If n' -■ n 
but s and s' different, we see writing out its value from (2) that the 
result is again zero. If 7i = n, s' = s, the result is 

'27r 



f 2t r2 

Jo Jo 



(1 - cosa - l3)"cos8acossPdadp. 



We may expand ( 1 - cosa - ^)" in cosines of multiples of a - /i, 
and the only term which contributes anything to the integral is that 
containing cos s(a - jS). 

Expanding by ordinary trigonometry, and integrating, we get 

(5) A( - i)v ^''- = ( - ly „_ ("i^y-M 

v/2"-i^ / (^_^,)|(^^^)| V \2n -{- l)(u - s)\(n + s)\ 

In order to obtain expansions of harmonics in polar co-ordinates, 

(ax + by-\- cz)" becomes l-^ — ^ ^ \^~V) ^cz\ 

= ^l^y(a^ + l7;-2t^)"say. 

Now taking a^ + — rj as a single term, we may expand this by 

a 

the Binomial Theorem, and pick out those terms in which the 

power indices of ^ and t] differ by s. 

Doing this we get (a^ + — rj - 2i«)" 

a 

8l(n - s)i 

Q^*^ being the function given in Thomson and Tait's Natural Philo- 
sophy, Vol. I., p. 205. Taking as a particular case the form (torcosa 
+ tysina -f zy, we easily find from this 



i tajcosa -\- y^ma -\- zYcossada = -— ■ —- r , ^JJ'cos5</> 

Jo ^ s\{n - sy 

Another mode of expansion leads to other well-known forms. 
Observe (a^ -|- — t) - 2iz)" may be written in either of the forms 
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V 

Expanding these in powers of a by Taylor's Theorem, we get the 
expressions for the harmonics in terms of differential co-efficients of 

The expansion of the Biaxal surface Harmonic may be deduced 
at once from (4) and (5). 

The definite integral expressions for the elementary harmonics 
are particular cases of expressions as the sum of a finite number of 
terms. 

Thus let V — (tarcosa + tt/sina + z)" = ^Hq + Hjcosa + H2COs2a + . . . . 

+ Kjsina + K2sin2a + . . . . 

and take a^, a^ a^ a series of angles in equi-different progression, 

the common difference being 27r/p. 

Then Vj = (Lxcosa^ + tysinaj + ;$;)"== J Hq + HjCOsaj + HaCOSOg + <fec., 
and it is easy to show that 

V,cos."<ai + V^cosa^ + V^cossa^ = J— H,, 

the other terms on the right hand side vanishing, provided 'p is not 
less than 2w+ 1. 

In this way we may express any harmonic as a sum of sectorial 
harmonics. 

In particular, 

Taking /? = oo , we get the definite integral 

— (tx'cosa + t2/sina + zfd(x 

as before. 

We may interpret this integral in a way that will lead us natu- 
rally to one or two others of similar form. 

Suppose we have a function of a:, y, 2;, the variables being con- 
nected by a relation, say 7? ^- 'f -^^ z^ = p^ \ so that to any particular 
value of z correspond in general an infinite system of values of a:, y, 

say a:i2/i, x^-i, &c. 

Then, obviously, the sum of all the values which the function 
takes for this particular value of z^ say /(a;,, 2/1, z)^rf(x^ 3/2, z)+ <tc. 
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(multiplied by a suitable infinitesmal) is a function of z and p alone. 
In particular cases, of course, it may be zero, or infinite. 

To treat the matter perfectly generally, we should have to con- 
sider cc, 2/, 2 as complex variables, but we may get a sufficient view 
under the following restrictions, viz. : we consider z as real, but x and 
y capable of taking values either purely real or purely imaginary. 

We may consider surface values, i.e., suppose a^ + y^ + z^ = l. 

Take then the harmonics (z + ta)**, (z + toj)"""^, both of which by 
(1) give us surface harmonics of the n^^ degree, and we get the fol- 
lowing cases : — 

I. z<l 



(i) x= Jl -z^cosa 
2/= ^1-Aina 

r2ir 

Harmonics are (^^ + 1 Jl - z^cosayda 

Jo 

•2t 



J 



^ {z + Lj\-z'cosa)—'^^ 



(ii) x= Jl -z^co&ha 

2/ = ^1 - 2;^tsinAa 

r+30 

Harmonic \^(^z + tj]- z^co^a)-^-^ da 



II. z>\ 



(i) x = cJz^-lcosa 
y=zi Jii?- Isina 
'2ir 



r2ir 

Harmonic (a- ^2;^ - lcosa)'»(/a 

Jo 



1 



2ir 

(z - Jz^ - lcosa)~"~^c?a 



(ii) x^ijz^-lco^a 

y = Jz^ - IsinAa 

r+00 

Harmonic \^{z- ^a* - lcosAa)-"-^cfa. 

These are well known expressions for Zonal Harmonics of the 
first and second kinds. By similar considerations we might also 
obtain integrals for the harmonics which are not zonal. 
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We go on to sketch briefly the application of the preceding 
methods to functions of p variables. 

The function v = {aiXi + a.p:.2+ ..0^0;^)" satisfies the equation 

The fundamental proposition about the integral of the product 
of two harmonics can be extended to the general case. To avoid 
circumlocution, we speak of points, lines, surfaces, volumes, by a 
well understood extension of the language of Analytical Geometry. 

Now we may extend Green's Theorem to the case of p variables, 
the proof following exactly the same lines as in the case of three 
variables. The only difficulty is about the meaning of the element 
of surface, but we may define it in this way. 

Let the surface be F(a;i, a;2...) = 0, and transform the variables 
linearly and orthogonally to a new set ^1, ^2> <^c-> of which ^1 = IjXi + 
l^i+ &c., where ^, l^, &c., are proportional to the values of 

— — , — — , &c., at the point Xi, x^', <fec., and are such that /i^ + 1^ + 
dxj 0X2 

• • • ^^ X • 

We have d^id^^ • • • d^^ = dx^dx^ . . . dx^y and we take d^^ • • • d^^ as 
the element of surface. We have then, putting dY for dx^dx^^.^dXp^ 
and dS for the element of surface 

We may, as before, apply this to finding the surface integral of 
the product of two harmonics (aiX^ + a^ +...)", {biX^ + 62^ +...)„ 
over the surface ajj' + ajj* + . . . = p\ 

The only difference in the general case is that the element of 
surface varies as r^^^ instead of as r^. Unless m = n, the integral 
vanishes. If m = w, the result is 

(remembering that -_ ( dxjdxo. ..dx„= Pf^ ^ Y 

^ dfjj r(?+i)' 



(6) P'"^'^''-^f^^^^{aA^aA..,Y. 



Let us now find an expression for the rational integral harmonic, 
which is a function of Xp and x^^ + cc^" + • • • + a^/. 
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I^ /i) ^21 • • 'fp-\ be real quantities such that the sum of their 

squares = 1 , then {f^x^ +/a + ■\'fp-iX^^ + ta-^)" is a harmonic, 

and if we take the sum of functions of this form for all real values 
of the f% such that the sum of their squares = 1, we shall get what 
we want. This mode of derivation of the symmetrical harmonic is 
somewhat different from our former method, but it is easily seen to 
lead to the same result. 

Consider first the integral 



w 



'"{fv^i^r +/_iiK^i + ta;^)"#i df^ 



i>-i 



with the integration extending to all real values of the/'s for which 
the sum of their squares is not greater than />-. 

Change the variables ^, f^ &c., by a linear and orthogonal sub- 
stitution to a new set, </>,, </>2, &c., of which 

^(x^ + . . . + a;Ai)^i =fv^\ + +./;>- i-V-i 

The integral is then 



11 



• • ( ^x{ + . . . a!^ii</>i + iXj^)"d(t>id<f>., . . . (/</>^i 



and the limits are not changed. 

Integrating with respect to </>.2 .•<^^>_i first, and disregarding 
numerical factors, we get 



J 



!f>° - </>i')^( W + -^A-^, + ue^rd<p,. 



Now differentiate this with respect to />, and in the result put 

This evidently gives us what we started to find. 
If we put </>i = cos^, the result takes the form 





Now the surface integral of the product of this and a harmonic 
(a^Xi + . . . + ctjjXpY is obtained by substituting a^ for jCj, a^ for x^, &c. 
in the above expression and multiplying by a numerical factor K. 

This integral is .-. K(ap)" r'^( 1 - cos^)"sin''-3^c?^ = M(ap)«, where 

Jo 

M is a numeric, easily found. 

By adding any number of harmonics of the typical form, we get 
a result which we may write 
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(7) [q„V„c?S = MV;, where V„' is the value of V„ at the pole 
of Q. 

The surface value of the above symmetrical harmonic is 



(x+ Jx" - loo&eYsin^-^Ode. 


We may consider this as a function of the single variable x, and 
it is easy from the above results to show that 



j: 



^}Q„Q«.(l-a^)^^aJ = 



when m, n are different, but = a certain finite numerical quantity, 
(jasily found, when m = n. 

The functions we have arrived at include, of course, as particular 
cases, simple and zonal Harmonic Functions. They are discussed at 
the end of Vol. I. of Heine's Ktigelfunctionen, 

The method of this paper is not mentioned in Heine, but I have 
found since I had it worked out that it is not new. In Professor 
Cayley's collected Works, Vol. I., page 397, will be found a short 
paper in which he proves (6) and (7). His proofs, however, are 
quite different from those I had arrived at, and have given above. 



A Method of Teaching Electrostatics in School. 
By J. T. Morrison, M.A., B.Sc. 

The object of the paper was to suggest for the teaching of electro- 
statics a leading idea, which should readily co-ordinate all the facts, 
introduce no misleading infer6Itces, and guide the course of learners 
in the direction of the most recent investigations — in all which 
respects the notion of attraction and repulsion is at least a partial 
failure. The leading idea or fact referred to is, that almost all 
electrostatic distributions, however complex, can be analysed into 
one or more repetitions of a certain simple system, which is called 
in the paper " an electrostatic system," and which may be described 
as follows : — Two equally and oppositely electrified conducting sur- 
faces, facing each other, separated by any dielectric, and insulated 
from each other. A complete study of one system of this kind, and 
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of the very simple ways in which the establishment of one such 
system often necessitates the establishment of others, is therefore 
the fundamental study of electrostatics. 

The following order may be adopted : — 

After getting pupils to discover that there are two kinds of 
electrification, and to find out the difference between conductors and 
non-conductors, the investigation of an " electrostatic system " is 
begun. 

A. Electric Quantity and Distribution in " an electrostatic " system. 

Chief Result : — The surfaces are equally and oppositely 
electrified, and to each small part of one surface may be shown 
to correspond a small equally but oppositely electrified part of 
the other surface. 

An important subsidiary principle may here be introduced 
— namely, the constancy of the quantity of electrification in a 
set of bodies that do not communicate with bodies that are 
outside the set. With the help of this principle, various dis- 
tributions should be analysed into their constituent " electro- 
static systems," and four important cases were so analysed in 
the paper. 

B. Mass-Motions in an " electrostatic system" 

Chief Result (a) ; — There is a mass-acceleration of each 
part of the two surfaces, the direction of the mass-acceleration 
being perpendicular to the surface. 

Definition of unit of electrification. 

(6) Mass acceleration on surface carrying unit positive 
electrification = - 4Tr(r. 

C. Energy and Potential Difference in an " electrostatic system,*^ 

Definition of Potential. Equipotential Surfaces. 
Chief Result : — Energy = J Quantity x Potential Differ- 
ence. 

Seat and Nature of Electrostatic Energy. 
Capacity, and Specific Inductive Capacity. 

D. Electric Discharge of a System, 

Electric Strength. 

Electric Oscillations. Transformations undergone by the 
Electric Energy. Electrical Radiation. 
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There is no room to indicate the experimental means by which 
the method may be executed, nor to show how the various instru- 
ments may be introduced, their laws discovered, and constants cal- 
culated ; but these matters will be fairly obvious. 



Eighth Meeting^ \^th June 1890. 



R. E. Allaedice, Esq., M.A., Vice-President, in the Chair, 



Note on the Orthomorphio Transformation of a circle into 

itself. 

From a letter by Professor Cayley. 

"The following is, of course, substantially well-known, but it 
strikes me as rather pretty : — to find the Orthomorphic transforma- 
tion of the circle 

ar* + y2-l=0 

into itself. Assume this to be 

1 -f- C(ac -I- ly) 
Then, writing A', B', C for the conjugates of A, B, C, we have 

AYaj - ly) + B' 

^ ^ \+C\x-iy)' 
and then 

o o ^ AAXx^ + f) + ABXa; + Ly) + A-B(a; - ty) + BE' 

'^' l + C{x + Ly) + C'{x-Ly) + CC\a^ + f) ' 

which should be an identity ior x^ + y^=l, a;? + y? = 1. 

Evidently = AB', whence C = A'B, and the equation then is 

1+AA'BB' = AA' + BB', 
i.e., (1-AA')(1-BB') = 0. 

But BB' = 1 gives the illusory result 

Xi + iyi = B, 
therefore 1 - A A' = 0, 

and the required solution thus is 

""'^'^''l + AB^x + LyY 
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where A is a unit- vector (say A = cos A, + tsinA.) and B, B' are conju- 
gate vectors. Or, writing B = 6 -H t/3, B' = 6 - t^, the constants are 
k, b, 13 ; 3 constants as it should be." 



Quaternion Synopsis of Hertz' View of the Blectrodynamioal 

Equations. 

By Professor Tait. 



Note on Menelaus's Theorem. 
By R. E. Allaedice, M.A. 

§ 1. The object of this note is, in the first place, to show that 
Menelaus's Theorem, regarding the segments into which the sides of 
a triangle are divided by any transversal, is a particular form of the 
condition, in trilinear co-ordinates, for the collinearity of three points; 
and, in the second place, to point out an analogue of Menelaus's 
Theorem in space of three dimensions. 

§ 2. In the usual system of areal co-ordinates, the a:-co-ordinate 
of P (fig. 52) is APBC/AABC, that is PD/AD. Now let D, E, F, 
be three points in BC, CA, AB, respectively, dividing these sides in 
the ratios /i/mj, ^/^n^, hh^ j then the co-ordinates of D, E, F, are 
proportional to (0, mj, /,), (Zg, 0, m^), {nt.j, Zj, 0). Hence the condition 
that D, E, F, lie on the straight line Ac + B// + Cs; = is 



m^ li 
L, m 
in.. L 



= 0, 



that is, ^1 /o ^3 + m^m.2m.i = 0, which is Menelaus's Theorem. 

§ 3. In space of three dimensions we may use the corresponding 
system of tetrahedral co-ordinates, and obtain a theorem analogous 
to that of Menelaus. 

Let BOD (fig. 53) be one of the faces of the tetrahedron ; and 
put a^ = PB'/BB' - APCD/BCD, Og = PCVCC = APDB/ACDB, etc. 
Then the co-ordinates of P, Q, R, S, points in the four faces of the 
tetrahedron, may be written (0, a.2, ag, a^), (6i, 0, 63, 64), etc. ; and 
the condition that these four points be coplanar is 
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a., a.^ a^ 


= 0, 


b, b., b. 




r'l c, C4 




d, il, (i, 





where «2> ^3> ^^^^ «4 "^^y ^^ taken to be the three areas into which 
the point P divides the face BCD ; and this condition is the analogue 
of Menelaus's Theorem for space of three dimensions. 



Historical notes on a geometrical problem and theorem. 

By J. S. Mackay. M.A., LL.D. 

The problem is 

Between two sides of a triangle to inflect a straight line which shall 
he eqiuil to each 0/ the segments 0/ the sides between it and the ba^e. 

This problem was brought before the Society at the January 
meeting in 1884, and a solution of it by Mr James Edward will be 
found in our Proceedings, Vol. II , pp. 5-6, a second by myself in 
Vol. II., p. 27 (10th April 1884), a third by Mr R. J. Dallas in 
Vol. III., pp. 41-2 (9th January 1885). Solutions of a slightly 
more general problem were also given by myself in Vol. III., pp. 
40-1, and reference made to the Educational Times, Vol. 37, p. 328 
(1st October 1884), and to Vuibert's Journal de Mathematiques 
itUmentaires, 9® ann^e, p. 45 (15th December 1884). 

I have since found that the more general problem was proposed 
by Monsieur J. Neuberg in the NouveVe Correspondance Mathd- 
inatique. Vol, I., p. 110 (1874-5), and solved by him in Vol. II, p. 
248 (1876); and quite recently I have discovered the first problem 
to go as far back as 1773-4. Here is how it occurs. 

In the Ladies^ Diary for 1773, Mr Thomas Moss proposes for 
solution the following : — 

TJie difference of the sides including a known angle of a plane 
triangle being given, and also the sum of one of those sides and that 
opposite the given angle, to construct tlie triangle* 

In 1774 the question is thus answered by Mr John Turner : — 

^^ Analysis. Suppose the thing done, and that ABC is the 



* Thomas Leybourn's Mathematical Questions proposed in the Ladies^ Diary, Vol. 
II., p, 377 (1817). 
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triangle, of which are given the angle A, the sum AC + CB, and the 
difference AB - AC. Then also will be given AB + BC, it being 
evidently equal to AC + CB + the given difference AB - AC ; and 
therefore if AD and AE be made respectively equal to the given 
sums AB + BO and AC + CB, the triangle ADE will also be given ; 
and then between the sides of this given triangle we have only to 
apply BC so as to cut off BD and CE each equal to it ; whence 
appears this easy 

" Construction. Form (Fig. 44) the given angle DAE and make 
DA and AE equal to the two given sums of each side and the base ; 
draw any line FG parallel to DE, and make EH and HI each equal 
to DF ; then draw EIB, and lastly BC parallel to HI cutting off 
the triangle ABC required." 

The demonstration is omitted as unnecessary. 

The theorem is 

The middle points of the three diagonals of a complete quad- 
rilateral are collinear. 

It was Camot (Essai sur la tMorie des transversales, 1806, § 6) 
who bestowed the name complete quadrilateral on a system of four 
straight lines no three of which are concurrent, but he does not seem 
to have been aware of this property of the three diagorals. In 
Cremona's Elements of Projective Geometry the property is ascribed 
to Gauss (Collected Works, Vol. IV., p. 391), and Dr Baltzer says, 
" Diese Eigenschaft des Vierecks ist von Gauss 1810 bemerkt worden 
(v. Zach monatl. Korresp. 22, p. 115)." 

In the Ladies^ Diary for 1795 the following question is proposed 
by Mr J. T. Connor, Lewes Academy, and in 1796 two solutions of 
it are given, the first by Major Henry Haldane of the Royal 
Engineers, and the second by Mr John Ryley of Leeds. 

In a triangle ABCy drawing any two lines AE, BD from, the 
extremities of one side, to terminate in the other tuH> sides; and 
tliereby fof^m a trapezium CDKE ; if the diagonals DE and CK of 
tJiat trapezium he bisected in the points M and N, and the first side 
of the triangle in the point P, these three points P, M, iV viill be in 
a straight line. Required the demonstration,^^ ^ 

* Thomas Leyboum's Mathematical Questions proposed in the Ladies^ Diary y Vol. ^ 
IIL, p. 291(1817). 
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